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The  structure  and  properties  of  Teflon-bonded  fuel  cell 
electrodes  and  their  components — platinum  black  and  Teflon 
30  particles--has  been  investigated  in  detail  and  compared 
with  previously  published  works.   Recent  experiments  em- 
ploying scanning  electron  microscopy  and  energy  dispersive 
X-ray  analysis  indicate  that  the  Teflon  particles  are  pro- 
late spheroidal  aggregates  of  approximately  spherical  Teflon 
molecules.   In  a  commonly  used  electrode  composition  these 
Teflon  aggregates  are  coated  with  one  or  more  layers  of 
catalyst  crystallites.   The  dimensions  of  the  Teflon  aggre- 
gates are  approximately  1500  x  3000  A,  the  diameter  of  the 
Teflon  molecules  ranges  from  200  to  about  300  A,  while  the 
diameter  of  the  catalyst  crystallites  appears  to  be  30  to 
40  A. 

Electrode  reactions,  ionic  hydration,  water  transport, 
and  their  effects  on  the  behavior  of  fuel  cells  have  been 
studied.   These  studies  suggest  that  hydration  of  potassium 
and  hydroxyl  ions,  which  has  been  ignored  in  previous  studies. 


XV 


causes  some  change  in  the  concentration  of  electrolyte  and 
its  gradient  in  the  electrolyte  matrix  and  the  electrode 
layer. 

On  the  basis  of  the  new  knowledge  of  the  electrode 
structure,  the  properties  of  its  components,  and  its  geo- 
metry a  mathematical  model  has  been  developed  to  enable  one 
to  predict  the  current  distribution  in  the  electrode  layer. 
These  calculations  show, as  one  might  expect,  that  at  high 
current  densities   the  great  majority  of  the  current  is 
generated  in  a  very  thin  layer  of  the  Teflon-bonded  catalyst 
immediately  adjacent  to  the  electrolyte  matrix.   For  an 

oxygen  electrode  operating  at  a  current  density  of  one  am- 

2 
pere  per  cm   more  than  half  of  the  current  is  generated  in 

the  20  %  of  the  catalyst-Teflon  layer  nearest  the  electro- 
lyte matrix.   At  low  current  densities  current  generation 
is  relatively  uniform  throughout  this  layer. 
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CHAPTER  I 
INTRODUCTION 

1.1   General  Description  of  Fuel  Cell  Structure 
and  the  Alkaline  Hydrogen-Oxygen  Fuel  Cell 

A  fuel  cell  is  a  direct  energy  conversion  device 
generating  electricity  by  oxidation  of  a  fuel  at  one 
electrode  and  reduction  of  an  oxidant  at  another  elec- 
trode.  Thus,  chemical  energy  is  converted  directly  into 
electrical  energy,  by  which  method  the  fundamental  limi- 
tation of  the  Carnot  cycle  is  avoided.   In  principle, 
therefore,  higher  efficiencies  can  be  attained  in  fuel 
cells  than  those  attained  in  heat  engines.   Hydrogen, 
hydrocarbons,  and  other  compounds  can  be  used  as  fuels 
while  oxygen  is  most  frequently  used  as  the  oxidant  either 
pure  or  as  air. 

Among  the  various  kinds  of  fuel  cells  which  have  been 
developed,  the  hydrogen-oxygen  system  in  alkaline  electro- 
lyte is  probably  the  most  advanced  and  widely  used.   How- 
ever, its  technology  has  far  outpaced  our  understanding  of 
its  behavior  and  operative  mechanisms. 

Since  gaseous  hydrogen  and  oxygen  are  used  as  fuel 
and  oxidant,  gas  diffusion  electrodes  are  used  for  both 
anode  and  cathode.   The  electrodes  are  usually  made  by 


depositing  a  mixture  of  catalyst  and  binder  on  a  wire 
screen  which  serves  both  as  an  electron  collector  and  a 
mechanical  support.   If  Teflon  is  used  as  the  binder,  it 
can  confer  a  degree  of  hydrophobicity  to  the  electrodes 
which  makes  it  possible  for  gaseous  reactants  to  diffuse 
readily  to  reaction  sites.   Between  the  anode  and  the 
cathode  there  is  an  electrolyte  matrix  which  provides 
electrical  insulation  between  the  electrodes  as  well  as 
a  porous  medium  through  which  diffusion  and  migration 
can  occur  under  the  influence  of  concentration  and  poten- 
tial gradients.   An  electrolyte  reservoir  is  also  used  to 
hold  extra  electrolyte  and  to  assist  in  evaporation  of 
water  produced  by  the  overall  cell  reaction;  it  is  usually 
a  sintered  metal  plate.   A  schematic  diagram  of  a  fuel  cell 
is  shown  in  Figure  1.1. 

The  electrochemical  reactions  occurring  in  this  fuel 
cell  are  usually  written  as  follows: 


Anode  :  2    E^      +      A    0H~         =   4  H2O   +   4  e 

Cathode  :         O2   +   2  HO   +   4  e   =   4  OH" 
Overall  Cell  Reaction  :    2  H2   +   O2   =   2  H2O 


However,  the  ions  are  hydrated  and  the  actual  reactions 
should  be  written  as  follows: 


Anode  :     2  H2   +   4 [OH  •SH2O]         =   4(l+s)  H2O  +  4  e 
Cathode  :   O2   +   2 (l+2s)  H2O   +   4  e   =   4[OH~-sH20] 
Overall  Cell  Reaction  :    2  H2   +   O2   =   2  H2O 

where  s  =  the  hydration  number  of  the  hydroxyl  ion. 


"2  ^" 
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Figure  1.1:   Schematic  Diagram  of  an  Alkaline  H^-O^  Fuel  Cell, 


since  water  is  produced  at  the  anode  and  consumed  at  the 
cathode,  both  in  the  formation  of  hydroxyl  ions  and  in 
their  hydration,  water  should  diffuse  from  the  anode  to 
the  cathode,  while  the  hydrated  hydroxyl  ions  should 
migrate  and  diffuse  from  the  cathode  to  the  anode.   The 
product  water  must  be  removed  from  the  cell  to  prevent 
electrolyte  dilution.   This  implies  more  complicated  trans- 
port processes  than  the  case  when  there  is  no  hydration  of 
hydroxyl  ions. 

1 . 2   Thermodynamics  and  Kinetics  of  Electrodes 

When  a  chemical  reaction  is  carried  out  electrochemically , 
the  work  done,  W,  is  given  by 

AH  =  AE  +  P  AV 

=   Q  -  W  +  P  AV  (1. 2-1) 

where  ah  =  enthalpy  change  of  the  system.  Joule/mole; 

^,E   -    internal  energy  change  of  the  system.  Joule/mole; 

P  =  pressure  of  the  system,  atm; 

3 

_1V  =  volume  change  of  the  system,  cm  ; 

Q  -    heat  absorbed  by  the  system.  Joule/mole. 

W  is  not  only  the  work  of  expansion  done  by  gases  produced, 
but  also  the  electrical  work  involved  in  transporting  elec- 
tric charge  around  the  external  circuit  from  the  anode  to 
the  cathode  at  which  the  potentials  are  V     and  V    ,  re- 
spectively .   The  maximum  electrical  work  that  can  be  done 


per  act  of  an  overall  reaction,  W ,  carried  out  in  a  cell 
and  which  involves  transfer  of  n  electrons  is 


W'   -  n  e  (V     -  V    )  (1.2-2) 

el  r ,  c     r ,  a 


for  a  hypothetical  case  in  which  the  internal  resistance  of 
the  cell  and  overpotential  losses  are  negligible.   Conver- 
sion to  a  molar  basis  by  multiplying  by  Avogadro's  number 
gives 

W  ,  =  n  F  (V     -  V    )  (1.2-3) 

el  r, c     r ,a 

where  F  =  Faraday's  constant,  equal  to  96,487  coulombs/ 

equivalent. 

Since  the  only  forms  of  work  involved  in  the  operation 

of  the  electrochemical  cell  are  electrical  work  and  work  of 

expansion, 

W  =  W  ,  +  P  AV  (1.  2-4) 

el 

If,  in  addition,  the  process  is  carried  out  reversibly, 

Q  =  T  AS  (1.2-5) 

where   AS  -    entropy  change  of  the  system,  Joule/mole-°K; 
T  =  temperature  of  the  system,  °K. 

Using  Equations  (1.2-1),  (1.2-3),  (1,2-4),  and  (1.2-5), 
it  follows  that 


AH  =  T  AS  -  n  F  (V     -  V    )  (1.2-6) 

r ,c     r ,a 


By  comparison  of  Equation  (1.2-6)  with 


AG  =  AH  -  T  AS  (1.2-7) 

where  AG  =  free  energy  change  of  the  system,  which  holds  for 
an  isothermal  process,  we  get 


AG  =  -n  F  (V     -  V    )  (1.2-8) 

r  ,c     r , c  ' 


Writing  V     -  V     =  E  ,  we  get 
r ,  c     r ,  a     r      ^ 


AG  -  -n  F  E^  (1.2-9) 


where  E   is  the  electromotive  force  of  the  reversible  cell. 
Thus,  an  ideal  fuel  cell  yields  electrical  energy  equivalent 
to  the  free  energy  change  of  the  reaction. 

In  the  normal  operation  of  a  fuel  cell,  the  electrical 
energy  must  be  supplied  at  a  significant  rate,  and  under 
these  conditions  reversibility  cannot  be  maintained.   In 
practice,  the  cell  potential,  E,  of  an  electrochemical  device 
working  at  an  appreciable  rate  (current)  is  given  by 


E  =  E   -  En  (1.2-10) 


where  En  is  the  sum  of  all  of  the  overpotentials  that  exist 
at  the  electrodes,  electrode-electrolyte  interfaces,  and 
ohmic  losses  in  the  electrolyte.   The  factors  that  affect 
the  working  voltage  of  a  cell  are  the  activation  overpoten- 
tial,  ohmic  losses,  and  mass  transfer  (concentration)  over- 
potential.   The  activation  overpotential  results  from  the 
slowness  of  one  or  more  of  the  reaction  steps  in  either  one 
or  both  of  the  electrode  reactions.   A  typical  plot  of 


terminal  cell  voltage  E  vs.  current  density,  i,  of  an 

operating  fuel  cell  is  shown  in  Figure  2.   The  open-circuit 

voltage  E    is  often  (usually  for  fuel  cells)  less  than  thermo- 
^    oc  ■' 

dynamically  reversible  potential  E^.  because  of  the  inter- 
ference caused  by  possible  reactions  due  to  impurities.   The 
effects  of  impurities  generally  decrease  as  the  current  den- 
sity increases.   It  is  always  found  that  there  is  a  sharp 
initial  decrease  of  the  working  potential  E  with  i,  as  shown 
by  section  AB  of  the  curve  in  Figure  1.2.   This  type  of 
behavior  is  characteristic  of  highly  irreversible  processes 
and  is  attributed  to  activation  overpotential. 

The  linear  portion  BC  in  Figure  2  corresponds  to  a 
relatively  high  current  density  region  in  which  the  decrease 
in  cell  potential  with  increase  of  current  density  is  due 
principally  to  ohmic  losses  in  the  electrolyte.   When  porous 
electrodes  are  used,  not  only  the  layer  of  electrolyte  in 
the  electrolyte  matrix  but  also  the  electrolyte  within  the 
pores  contributes  to  the  electrolyte  resistance. 

At  sufficient  high  current  densities  (i.e.,  high  local 
reaction  rates) ,  most  heterogeneous  reactions  pass  into  a 
region  where  the  rate  is  controlled  by  the  rate  of  transport 
of  reactants  to,  or  products  from,  the  electrodes.   Region 
CD  in  Figure  1.2  illustrates  this  effect,  which  is  known  as 
the  concentration  (or  mass  transfer)  overpotential. 

In  the  actual  operation  of  fuel  cell,  the  current  I, 
after  considering  the  overpotentials ,  is  given  by 


V 


^'C   ,, 


V 


r,a 


iim,c        lim,a 


lim 


Current  Density,  i 


Figure  1.2:   A  Typical  Cell  Potential-Current  Density 
Relation  for  a  Cell. 


I  =    -^ 1 2hSl  (1.2-11) 

R 
e 


where   n  ,   -    ohmic  voltage  loss,  volt; 
ohm  ^ 

R     =  resistance  of  the  external  load,  ohm; 
e 

V^;  and  Va=  electrode  potentials  at  the  cathode  and 

the  anode,  respectively,  which  are  given  by 


V   =  V     -    T-r] 
c    r,c     c 

=  V    -  n    4.  -  n  (1.2-12) 

r,c     c,act     c,conc 


V  ==  V    +  Zn 

a    r,a     a 

=  V    +  n    ^  t  n  (1.2-13) 

r,a     a, act    a, cone 


where  the  suffixes  "act"  and  "cone"  denote  activation  and 
concentration  overpotentials ,  respectively. 

In  many  electrochemical  reactions,  the  current  at  a 
particular  electrode  is  found  to  vary  exponentially  with 
potential  across  the  metal-solution  interface,  and  the  re- 
lationship between  potential  and  the  rate  constant  k  of  an 
electrochemical  reaction  can  be  expressed  empirically  by 

k  =  k°  exp  (-3VF/RT)  (1.2-14) 

where    k      =    rate    constant    in    the    absence    of    a  potential  difference 
B      =   a    symmetry    factor,    0    <    6    5    1 
R      =    gas    constant,    8.314    Joule/gmole-°K 

The  reaction  velocity  in  electrochemical  kinetics  is 
related  to  the  current  density,  i,  which  is  the  current  per 
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unit  area  of  electrode  surface.   The  relation  between  re- 
action velocity,  v,  and  current  density,  i,  depends  on  the 
number  of  electrons  transferred  in  one  act  of  the  overall 
reaction,  n. 

i  =  n  F  V  (1.2-15) 

The  rates  of  electrochemical  reactions  are  a  measure  of 
the  net  reaction  rates  and  the  net  current  density  is  given 
by 

i  =  1  -  1  (1.2-16) 

where  i  is  the  forward  current  and  i  is  the  reverse  current. 
Therefore,  in  the  case  of  the  electrode  reaction  O  +  e  ->  R, 
the  net  current  density  may  be  expressed  by  (since  n  =  1) 


i  =  F 


kf  Cq  exp 


evF 


RT 


-  k   C^  exp 
r   R    "^^ 


(1-6) 


VF 
RT 


(1.2-17) 


where  'kr:    and  k   are  the  rate  constants  of  the  forward  and 
f       r 

reverse  reactions  in  the  absence  of  a  metal-solution  poten- 
tial difference  and  C„  and  C^  are  the  concentrations  of  the 

oxidized  and  reduced  species,  respectively.   The  exchange 

2 
current  density  is  related  to  the  velocity,  in  amp/cm  ,  of 

the  forward  or  backward  reaction  at  equilibrium,  in  which 

condition  the  net  current  density  is  zero  and  the  exchange 

current  density  i   is  given  by 


F  k°  Cq  exp 


3VF 


RT 


F  k°  C„  exp 

r   R    '^ 


<i-6)  if 


(1.2-18) 
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Using  Equation  (1.2-18)  and  the  overpotential ,  which 
is  equal  to  the  potential  difference  V  -  V  ,  Equation 
(1.2-17)  can  be  expressed  as 


1=1 


exp 


BnF 


RT 


-  exp  (1-3) 


nF 
RT 


(1.2-19) 


In  the  case  of  a  multielectron  transfer  reaction,  the  general 
form  of  the  current  density-overpotential  relation  is 


1=1 


exp 


a  riF 
a 

RT 


exp 


a^nF 


RT 


(1.2-20) 


where  a   and  a   are  the  transfer  coefficients  of  the  reactions 
a      c 

of  the  anode  and  the  cathode,  respectively,  and  depend  on 
the  mechanism  of  the  overall  reaction.   In  the  study  of 
mathematical  formulation  and  solution  of  electrode  models, 
this  equation  is  most  important. 

1.3   Advantages  of  Porous  Gas  Diffusion  Electrodes 


For  an  electrochemical  energy  conversion  system  to  be 
useful  for  practical  applications  it  should  be  able  to  pro- 
vide good  power  density.   Because  of  the  need  to  minimize 
the  activation  overpotential  at  desired  currents,  to  facil- 
itate mass  transfer  to  the  reaction  sites,  and  because  both 
the  fuel  and  the  oxidant  in  hydrogen-oxygen  fuel  cells  are 
gaseous,  porous  gas  diffusion  electrodes  are  used.   The  much 
greater  power  densities  afforded  by  porous  electrodes  as  com- 
pared to  planar  electrodes  is  due  not  only  to  their  high 
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rates  of  real  to  apparent  area  (i.e.,  much  larger  concentra- 
tions of  reaction  sites  are  available) ,  but  also  to  the  fact 
that  they  permit  much  higher  limiting  currents.   Another  ad- 
vantage of  porous  electrodes  is  that  the  catalyst,  in  the 
form  of  very  fine  particles,  can  be  dispersed  within  the 
pores,  thereby  reducing  considerably  the  quantity  of  catalyst 
used . 

The  difficulty  of  understanding  the  behavior  and  mecha- 
nisms of  electrode  reactions  in  porous  electrodes  arises 
mainly  from  their  extremely  complex  physical  structure  and 
the  consequent  difficulty  of  describing  them  geometrically 
so  that  mathematical  analysis  can  be  carried  out.   One  of 
the  principal  purposes  of  this  study  is  to  examine  the 
detailed  structure  of  Teflon-bonded  fuel  cell  electrodes  and 
to  relate  this  structure  to  their  performance  and  transport 
processes,  etc.,  so  that  we  can  use  the  information  in  under- 
standing and  evaluating  the  behavior  and  mechanisms  of  per- 
forming electrodes  and  in  designing  electrodes  capable  of 
improved  performance. 


CHAPTER  II 
THE  STRUCTURE  OF  TEFLON-BONDED  ELECTRODES 

2  . 1   Characterization  of  Electrode  Material 

Teflon-bonded  fuel  cell  electrodes  are  made  of  three 
different  materials;  catalyst.  Teflon,  and  fine-meshed 
metal  screen.   In  many  cases  platinum  black  is  used  as  the 
catalyst  because  of  its  excellent  catalytic  properties, 
i.e.,  its  high  exchange  current  density  in  alkaline  elec- 
trolyte.  It  is  a  very  fine  black  powder;  the  crystallite 
size  is  so  small  that  individual  crystallites  have  not  yet 
been  clearly  resolved  by  scanning  electron  microscopy,  i.e., 
the  crystallite  diameter  is  less  than  10  nm  (=100  A) .   How- 
ever, by  surface  area  measurements  using  the  BET  gas  adsorp- 
tion method  [1],  the  average  crystallite  diameter  of  a 
widely  used  platinum  black  (Englehard  fuel  cell  grade)  is 
calculated  to  be  about  40  ±  10  A.   In  its  as-received  state, 
the  catalyst  exists  as  agglomerates  of  1  micron  in  diameter. 

Teflon  30,  an  aqueous  emulsion  (DuPont) ,  is  usually 
used  for  the  Teflon  component  of  electrodes.   It  is  a  nega- 
tively charged  hydrophobic  colloid  containing  polytetra- 
f luoroethylene  (PTFE)  resin  particles,  0.05  to  0.5  micron 
in  diameter,  suspended  in  water  and  stabilized  with  Triton 
X-100,  a  nonionic  surfactant.   Teflon  is  very  hydrophobic 
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and  is  relatively  stable  up  to  temperatures  of  about  350°C. 
Electron  microscopic  examination  of  well-dispersed  pure 
Teflon  particles  reveals  that  although  the  size  varies  from 
500  to  5000  A  in  diameter,  the  majority  of  particles  are 
prolate  spheroid  in  shape,  the  minor  and  major  diameters 
being  about  1000  to  1500  A  and  2000  to  3000  A,  respectively. 
It  is  also  noted  that  a  single  Teflon  particle  appears  to 
be  an  aggregate  of  small  spherical  particles  (200  to  300  A 
in  diameter)  which  may  be  Teflon  molecules.   Teflon  30  is 
reported  by  the  manufacturer  to  have  a  mean  molecular  weight 
of  5  MM,  and  calculation  suggests  that  a  molecule  of  this 
molecular  weight  should  have  a  diameter  of  about  200  A. 
Several  scanning  electron  micrographs  of  well-dispersed  Teflon 
particles  are  shown  later  on  pages  28-29. 

The  fine-meshed  metal  screen  is  used  to  provide  mechan- 
ical support  for  the  catalyst  layer  and  to  serve  as  an  elec- 
tron collector  for  the  electrons  generated.   Gold-plated 
nickel  screen  about  100  mesh  is  used  for  most  alkaline  fuel 
cells. 

2 . 2   Surface  Area  Measurements 

2.2.1   Method  of  Measurements 

The  surface  areas  of  platinum  catalyst,  pure  Teflon, 
and  electrodes  were  measured  by  BET  method.   This  method 
involves  measuring  the  amount  of  gas  adsorbed  on  a  solid 
surface  at  a  temperature  close  to  the  boiling  point  of  the 
gas.   A  quantity  of  sample  is  weighed  into  a  sample  cell 
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and  pretreated  for  about  15  hours  at  200-220°C  while  a 
stream  of  pure  helium  flows  through  the  cell  to  remove 
moisture  and  adsorbed  species.   After  the  measurement  the 
sample  is  reweighed  to  determine  the  dry  weight. 

In  the  continuous  flow  method  developed  by  Nelsen 
and  Eggertsen  [2]  nitrogen  is  used  as  the  adsorbate  and 
helium  is  used  as  the  carrier  gas.   A  mixture  of  nitrogen 
and  helium  of  known  composition  is  passed  through  the 
sample  and  the  effluent  composition  is  monitored  by  thermal 
conductivity  measurements.   When  the  sample  is  cooled  in 
liquid  nitrogen,  the  adsorption  of  nitrogen  is  indicated  by 
a  peak  on  a  recorder  chart.   After  adsorption  equilibrium 
is  established,  the  recorder  pen  returns  to  its  original 
baseline.   The  sample  tube  is  allowed  to  warm  by  removing 
the  liquid  nitrogen  coolant,  causing  desorption  of  nitrogen 
and  producing  a  peak  on  the  chart  which  is  in  the  reverse 
direction  of  the  adsorption  peak. 

The  area  under  the  desorption  peak  is  a  measure  of 
the  nitrogen  adsorbed.   This  peak  area  was  calibrated  by 
injecting  a  known  amount  of  nitrogen  into  the  nitrogen-helium 
stream  to  give  a  peak  of  similar  magnitude  to  the  area  of 
the  desorption  peak.   The  amount  of  nitrogen  adsorbed  is  cal- 
culated by  comparing  these  two  peak  areas.   Adsorption,  de- 
sorption, and  calibration  are  done  at  several  relative  nitrogen 
pressure  ratios  in  the  range  from  0.05  to  0.3. 

After  computing  the  amount  of  nitrogen  adsorbed  and  the 
corresponding  relative  pressures,  the  amount  of  nitrogen  ad- 
sorbed in  a  monolayer,  V  ,  is  determined  from  the  BET  equation: 
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"-    '       '      *J^  (2.2-1) 


V,(P-P)    VC    P     vc 
ads    o         mom 


where  P  -    partial  pressure  of  the  adsorption  gas,  mmHg; 

P  -    saturation  pressure  of  the  adsorption  gas  over 
the  solid  sample  at  the  temperature  of  coolant, 
mmHg; 

V  ,   =  total  volume  (STP)  of  adsorbed  gas  on  the  sur- 
face of  the  adsorbent,  ml ; 

C   =  constant  expressing  the  net  adsorption  energy; 

V   =  volume  (STP)  of  adsorbed  gas  when  the  entire 
m 

adsorbent  surface  is  covered  with  a  monomolecular 

layer,  ml. 


When  the  left  hand  side  of  the  equation  is  plotted  as  the 
ordinate  versus  the  relative  pressure,  P/P  ,  a  straight  line 
with  slope  (C  -  1)/V   C  and  intercept  1/V   C  is  obtained. 
If  we  note  that 


=  V   ,    (2.2-2) 


Slope  +  Intercept    C  -  1     1      m 

V   C     V   C 
m        m 


then  we  can  determine  V   from 

m 


V   =  -, ^n—  (2.2-3) 

m    k,  +  k^ 


where  k,  is  the  slope  and  k„  is  the  intercept  of  the  straight 
line.   This  value,  multiplied  by  the  proper  factor  for  the 
area  covered  per  unit  amount  of  nitrogen,  gives  the  surface 
area,  and  when  divided  by  the  dry  weight  of  the  sample,  it 
yields  the  specific  surface  area. 
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2.2.2   Experimental  Results  and  Discussion 

Weighed  samples  of  platinum  catalyst  (Englehard  fuel 
cell  grade  platinum  black)  were  heated  in  a  stream  of 
flowing  helium  for  about  15  hours  at  temperatures  ranging 
from  100  to  400  °C.   Then  the  surface  area  was  measured  as 
described  in  Section  2.2.1.   Since  it  is  standard  procedure 
to  sinter  fuel  cell  electrodes  fabricated  with  Teflon- 
platinum  catalyst  mixtures  at  a  temperature  between  310  and 
320  °C,  the  surface  area  in  the  temperature  range  300-350  °C 
was  studied  carefully.   The  data  are  shown  in  Table  2.1  and 
in  Figure  2 . 1 

Table  2.1 

Dependence  of  Surface  Area  of  Platinum  Black  on 
Pretreatment  Temperature. 


Temp,  of  Pretreatment,  °C 


Surface  Area,  m  /g 


100 
150 
200 
250 
300 
310 
320 
330 
350 
400 


45.3 
46.0 
48.9 
49.6 
48.1 
44.2 
37.  8 
12.1 
5.25 
5.01 
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50 


40 


30 


20 


10 


100 


200  300 

Temperature    of    Pretreatment ,     °C 


400 


Figure    2.1; 


Dependence    of    Surface   Area   of   Pt-black   on 
Pretreatment   Temperature. 
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It  is  clear  that  a  large  decrease  in  surface  area 
occurred  when  the  pretreating  temperature  was  between  300 
and  350  °C,   There  is  a  gradual  increase  in  the  surface 
area  as  the  pretreating  temperature  is  increased  to  200  °C, 
but  the  slight  changes  in  surface  area  between  200  and  300  °C 
are  probably  not  significant.   However,  the  reduction  of 
surface  area  seems  to  have  already  begun  at  300  °C,  and  the 
advantages  of  sintering  the  Teflon  at  310  -  320  °C  must  be 
weighed  against  the  loss  in  surface  area.   On  the  other  hand 
the  surface  area  lost  on  sintering  may  be  lost  fairly  early 
during  the  life  of  an  electrode,  and  it  may  not  be  signif- 
icant for  the  entire  life  of  a  fuel  cell. 

2.2.3   Estimation  of  Total  and  Free  Surface  Areas  of 
Electrode  Components 

The  surface  areas  of  platinum  black  catalysts  and  Teflon 
particles  can  be  estimated  by  calculation  from  the  informa- 
tion obtained  on  their  sizes.   Assuming  that  both  catalyst 
crystallites  and  Teflon  molecules  are  spherical  and  their 
surfaces  smooth,  we  can  calculate  the  surface  areas  from 


S   =  (Surface  Area  of  a  Single  Particle)  X  (Number  of 
Particles  in  1  gram) 


TT  d^  X 


T  /'T  j3 
1/g  d  p 


whe 


-  6/pd  (2.2-4) 

2 

re   S   =  total  geometrical  surface  area  per  one  gram,  cm  /g; 

3 

p   =  density  of  the  component,  g/cm  ; 

d   =  diameter  of  the  component,  cm. 
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In  comparing  the  calculated  surface  areas  with  those 
measured  by  the  BET  method  the  exclusion  of  surface  area 
due  to  the  presence  of  neighboring  particles  should  be 
considered.   The  fraction  of  excluded  surface  area  can  be 
estimated  from  the  following  [3]: 

n  d 

^  =  2  (d   ;  d)  (2.2-5) 

o 

where   f  =  fraction  of  excluded  surface  area; 

n  =  coordination  number,  i.e.,  number  of  contacts 
per  particle; 

d   =  diameter  of  adsorbate  gas  used  in  the  BET  method, 
o  ^ 

cm; 

d  =  diameter  of  particle,  cm. 

Therefore,  by  using  Equations  (2.2-4)  and  (2.2-5)  and  noting 
that  the  specific  gravity  of  Teflon  is  2.22  and  that  of 
platinum  is  21.45  while  letting  the  molecular  diameter  of 
Teflon  and  Pt-black  range  from  200  to  300  A  and  10  to  60  A, 
respectively,  d   =  4.5  A  for  nitrogen,  and  n  =  4,  6  and  8, 
we  get  the  results  shown  in  Table  2.2. 

Comparing  the  results  with  the  surface  area  of  platinum 

2 

black  experimentally  measured,  45  m  /g,  we  can  estimate  that, 

assuming  the  packing  condition  of  platinum  black  to  be  random 
to  random-close,  the  average  size  of  a  single  catalyst  crys- 
tallite ranges  from  45  to  50  A  in  diameter.   Considering  the 
assumption  that  the  crystallites  are  uniform  smooth  spheres, 
we  can  see  that  the  actual  size  of  a  crystallite  can  be 
smaller  than  what  has  been  estimated,  i.e.,  ranging  from  30 
to  4  0  A  in  diameter. 
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Table  2.2 
Total  and  Free  Surface  Areas  of  Teflon  and  Pt-black. 

2 
Diaireter,   Calculated  Geometric   Free  Surface  Areas,  m  /g 

A       Surface  Area,  ra^/g    n=4   n=6    n=8 

200  135  129     126    123 

Teflon        250  108  104     102    100 

300  90  87      86     85 

10  279  106      19 

20  140  89      63     37.1 

Platinum  (25)43.7 

Black       30  93  69      57     44.5 

(35)43.5 
40  70  56      49     41.7 

50  56  47      42     37.5 

60  47  40      37     33.9 


2.2.4   Estimation  of  Catalyst  Layer  Thickness  on  Teflon 
Aggregates 

In  order  to  be  able  to  estimate  the  thickness  of  Pt- 
black  layer  on  Teflon  aggregates  it  is  necessary  to  know  the 
surface  area  of  Teflon  aggregates  and  the  packing  condition 
of  Pt-black  crystallites  on  the  surface  of  Teflon  aggregates. 
From  the  information  on  the  size  and  the  shape  of  Teflon  ag- 
gregates we  can  calculate  the  surface  area  of  Teflon  aggregates 
as  follows: 
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Surface  Area  of  1  g  of  Teflon  =  (Surface  Area  of  a  Teflon 
Aggregate)  X  (Number  of  Aggregates  in  1  g  of  Teflon) 

1 


TTr  , -,   ,   A  B   .    , 


X 


-6-  Vb  P  (1  -  ^). 


(2.2-6) 


where   d   =  major  diameter  of  a  Teflon  aggregate,  cm; 

d„  =  minor  diameter  of  a  Teflon  aggregate,  cm; 

/~~2  2 

e      =    eccentricity,    e    =  /  d^    -    d^    /   d^    ; 

r  -    roughness  factor; 

e   =  porosity. 


A 


B 


A 


Therefore,  substituting  d^  =  2500  A,  d„  =  1250  A,  r  =  1.5, 

A  n 

e    =    0.36  (for  random-close  packing),  and  p  =  2.22,  we  get 

2 

the  external  surface  area  of  Teflon  aggregates  as  43.3  m  /g, 

The  thickness  of  Pt-black  catalyst  layer  in  the  mix- 
tures of  Pt-black  and  Teflon  can  be  calculated  as  follows: 


Thickness  =  (Bulk  Volume  of  Pt-black  in  1  g  of  Mixture)/ 
(Teflon  Aggregate  Surface  Area  in  1  g  of  Mixture) 
1-x 


p(l-£) 


/  (Surface  Area  of  1  g  of  Teflon  Aggregates  x  x) 

(2.2-7) 


where   x  =  fraction  of  Teflon  in  a  mixture; 

3 
p  =  density  of  Pt-black  crystallites,  g/cm  ; 

e  =  porosity  of  Pt-black  catalyst  layer. 

Therefore,  with  an  assumption  that  the  packing  of  Pt-black 
crystallites  on  Teflon  aggregates'  surface  is  random,  which 
has  the  porosity  of  0.36  to  0.40,  we  can  estimate  that  an 
80:20  Pt-black-Teflon  mixture  has  a  catalyst  layer  thickness 
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of  70  A  and  a  50:50  mixture  has  a  catalyst  layer  thickness 
of  17  A.   Considering  the  fact  that  the  crystallite  size 
of  Pt-black  has  been  estimated  to  be  30  to  40  A  in  diameter, 
we  can  see  that  in  the  50:50  mixture  a  fairly  significant 
fraction  of  the  surface  of  Teflon  aggregates  is  not  coated 
with  Pt-black  catalyst,  although  there  can  still  be  continuous 
catalyst  connection  due  to  the  tendency  of  catalyst  particles 
collecting  in  the  crevices  first. 

2.2.5   Comparison  of  the  Calculated  and  Measured  Surface  Areas 
of  80:20  and  50:50  Pt-black-Tef Ion  Layers  and  Pure 
Teflon  Layer 

If  we  assume  that  the  surface  areas  of  the  components 

of  Pt-black-Tef Ion  layers  do  not  change  during  the  formation 

of  electrodes,  we  can  estimate  the  surface  areas  of  Pt-black- 

Teflon  mixture  layers  of  different  ratios  from  calculations 

by  using  pure  component  surface  areas.   Comparison  of  these 

calculated  surface  areas  with  the  measured  ones  is  shown  in 

Table  2.3. 

Table  2.3 

Calculated  and  Measured  Surface  Areas  of  Pt-black- 
Teflon  Mixtures  and  Pure  Teflon. 


Pt-black:Teflon  Calculated  Measured 

wt.  %  ratio        Surface  Area,  m  /g      Surface  Area,  m  /g 


80:20  57     (38)  22.0    ±    0.5 

50:50  74     (27)  19. 0+0.5 

Pure    Teflon  103  9.5    +    0.3 
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As  we  can  see  in  the  results  shown  above,  there  is  a 
significant  difference  between  the  calculated  and  the  mea- 
sured surface  areas  of  pure  Teflon.   Since  the  pure  Teflon 
sample  was  prepared  by  filtration,  pressing,  and  sintering 
for  the  measurement  of  its  surface  area,  there  must  have 
occurred  a  considerable  reduction  of  its  surface  area  during 
these  processing,  particularly  due  to  the  sintering  at  310°C 
which  would  cause  softening  and  binding  of  Teflon  aggregates. 

If  we  use  the  measured  surface  area  of  Teflon, 

2 
9.5  +  0.3  m  /g,  in  calculating  the  surface  area  of  the  mix- 
tures, we  obtain  the  results  shown  in  the  parenthesis  in 
Table  2.3.   Comparing  these  results  with  the  measured  surface 
areas  of  the  Pt-black-Tef Ion  layers  of  different  ratios,  we 
can  still  notice  considerable  differences  between  the  cal- 
culated and  measured  surface  areas.   This  can  also  be  attri- 
buted to  the  softening  and  bonding  of  Teflon  aggregates 
during  sintering,  which  can  cause  some  portion  of  Pt-black 
to  be  buried  in  them.   This  can  be  shown  by  the  fact  that  the 
calculated  surface  area  of  Pt-black  from  the  measured  surface 

areas  of  80:20  and  50:50  mixtures  and  pure  Teflon  is  about 

2 
27  m  /g,  which  would  mean  that  there  was  a  loss  about  40  % 

in  the  surface  area  of  Pt-black. 

2 . 3   Formation  of  Electrodes 

Fuel  cell  electrodes  can  be  made  in  several  ways;  one 
method,  based  on  information  derived  from  reports  and  tech- 
nical personnel  at  NASA  Laboratories  is  as  follows: 
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Platinum  black  is  weighted  out  and  dispersed  for 
several  minutes  in  water  containing  a  wetting  agent, 
Triton  X-100,  in  an  amount  equal  to  6%  of  the  weight  of 
catalyst  using  an  ultrasonic  agitator.   The  requisite 
amount  of  Teflon  30  emulsion  is  added  to  Pt-black  disper- 
sion and  the  ultrasonic  agitation  is  repeated  for  an  addi- 
tional several  minutes.   The  mixture  is  filtered  on  a  very 
fine  pore  filter.   After  the  catalyst  mixture  is  partially 
dried  at  room  temperature,  a  fine-meshed  metal  screen  is 
placed  on  the  catalyst  layer,  followed  by  laying  a  filter 
paper  on  the  screen  to  absorb  water  and  the  "sandwich"  is 
pressed  at  about  5000  psi.   The  filter  papers  are  then  peeled 
off,  and  the  electrode  is  sintered  in  air  at  310  °C  for  about 
20  minutes. 

2 . 4   Structure  of  Electrodes 

In  understanding  various  phenomena  occurring  in  an 
operating  fuel  cell,  it  is  essential  to  know  the  structure 
of  the  catalyst  layer  as  this  has  a  great  influence  on  the 
performance  of  the  cell.   Since  the  catalyst  layer  of  an 
electrode  consists  of  a  mixture  of  very  small  hydrophylic 
Pt-black  particles  and  hydrophobic  Teflon  particles,  their 
arrangement  in  the  layer  will  have  a  great  influence  on  the 
wettability  of  electrodes,  thereby  influencing  the  transport 
processes  within  electrodes.   The  relative  amounts  of  elec- 
trode components,  as  well  as  their  arrangement,  will  have 
influence  on  these  characteristics. 
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Most  of  the  published  work  on  the  structure  of  Teflon- 
bonded  fuel  cell  electrodes  [4,5]  is  based  on  the  view  that 
the  basic  structure  of  the  catalyst  layer  is  one  of  fairly 
large  catalyst  aggregates  more  or  less  surrounded  by  Teflon- 
particles,  the  two  materials  existing  as  interpenetrating 
lattices  which  provide  continuous  transport  paths  for  both 
liquid  and  gas.   This  view  is  illustrated  in  the  schematic 
diagram  below. 


TEFLON 
PARTICLES 


LIQUID  PHASE 
BRIDGES 


CATALYST 
AGGREGATE 


GAS 
SPACE 


Our  recent  experiments  employing  scanning  electron 
microscopy  and  energy  dispersive  X-ray  analysis  have  raised 
serious  questions  as  to  the  accuracy  of  the  view  described 
above.   Indeed,  they  suggest  very  strongly  that  the  struc- 
ture of  the  catalyst-Teflon  layer  is  opposite  to  it. 
Scanning  electron  micrograph  of  electrodes  and  of  disper- 
sions of  Pt-black,  Teflon,  and  their  mixtures  in  different 
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ratios  have  been  made,  and  these  reveal  several  significant 
details  of  the  structure  of  the  catalyst-Teflon  layer  as 
summarized  below; 

1)   The  Size,  Shape,  and  Substructure  of  Teflon 
Particles 

Scanning  electron  micrographs  (SEM)  of  pure 
Teflon  settled  on  a  glass  coverslide  from  a  dilute 
dispersion  revealed  that  Teflon  particles  are 
roughly  prolate  spheroidal  in  shape,  but  there 
is  some  evidence  that  these  particles  are  aggre- 
gates of  closely  packed  spheres  as  can  be  seen 
in  Figure  2.2.   The  substructure  of  the  Teflon 
particle  is  substantially  more  obvious  in  the  SEM 
of  catalyst-Teflon  mixtures  as  seen  in  Figures 
2.3  and  2.4,  which  are  scanning  electron  micro- 
graphs of  Pt-black-Tef Ion  suspensions  settled  on 
filter  paper  and  in  Figure  2.5  which  is  a  scanning 
electron  micrograph  of  a  fuel  cell  electrode.   The 
increased  visibility  of  the  substructure  may  be 
due  to  the  contrast  provided  by  catalyst  crystal- 
lites packing  in  the  crevices  between  Teflon  spheres. 

If  one  assumes  that  the  individual  molecules 
are  spherical,  calculations  by  using  mean  molecular 
weight  of  5  MM  and  density  of  2,22  show  that  the 
size  of  a  Teflon  molecule  ranges  from  200  to  300  A 
in  diameter. 


Figure  2.2:   Scanning  Electron  Micrograph   of  Pure  Teflon 
Particles.   Magnification:   50,000. 


Figure  2.3: 


Scanning  Electron  Micrograph   of  Pt-black- 
Teflon  Mixtures.   Magnification:   36,000. 
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Figure  2 . 4 


Scanning  Electron  Micrograph   of  Pt-black- 
Teflon  Mixtures.   Magnification:   72,000. 


Figure  2.5 


Scanning  Electron  Micrograph   of  a  Fuel  Cell 
Electrode.   Magnification:   30,000. 
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2)   Distribution  of  Pt-black  in  the  Catalyst-Teflon 
Layer 

It  has  been  noted  earlier  that  scanning  elec- 
tron microscopy  has  not  resolved  individual  catalyst 
crystallites,  but  electron  micrographs  of  catalyst- 
Teflon  mixtures  suggest  very  strongly  that  the 
catalyst  particles  coat  the  surface  of  the  Teflon 
aggregates. 

At  large  catalyst-Teflon  ratios   (Pt  >  80  wt.  %) 
multiple  layers  of  catalyst  particles  appear  to 
coat  the  surface  of  the  Teflon  aggregates  and  to 
fill  the  interstices  between  them.   This  phenomenon 
appears  to  be  further  substantiated  by  the  fact  that 
Teflon  particles  are  negatively  charged  and  Pt-black 
particles  are  positively  charged.   They  would, 
therefore,  be  attracted  to  each  other.   This  view 
of  the  catalyst-Teflon  layer  is  illustrated  in  the 
diagrams  below. 
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3)   Particle  Packing  and  Pore  Structure 

The  packing  arrangement  seems  to  be  very 
close  among  Teflon  molecules  in  an  aggregate,  and 
there  were  obviously  macropores  existing  in  a  more 
or  less  random  manner. 

In  order  for  porous  gas  diffusion  electrodes 
to  be  effective  there  must  be  an  appropriate  amount 
of  catalyst  which,  after  being  wetted  by  electro- 
lyte, provides  reaction  sites  for  the  heterogeneous 
electrochemical  reactions.   Moreover,  the  catalyst 
layer  must  also  provide  electron  conduction  paths, 
liquid  phase  mass  transfer  paths  for  reactants  and 
products,  and  gaseous  flow  paths.   When  the  propor- 
tion of  catalyst  in  the  catalyst-Teflon  layer  is 
too  great,  the  catalyst  will  form  a  more  or  less 
continuous  matrix  so  that  Teflon  particles  will 
exist  as  isolated  "islands,"  in  which  condition 
there  would  be  no  continuous  gaseous  channels  with- 
in the  catalyst  Teflon  layer.   This  would  cause 
very  limited  utilization  of  the  electrodes  owing 
to  the  large  gas  diffusion  resistance,  i.e.,  almost 
all  reaction  would  occur  near  the  electrode-gas 
interface,  resulting  in  a  sharp  decrease  in  elec- 
trode efficiency.   On  the  other  hand,  when  there 
exists  too  high  a  proportion  of  Teflon,  the  dis- 
tribution of  catalyst,  as  well  as  electrolyte  paths, 
may  be  discontinuous,  in  which  circumstance  the 
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reactions  would  occur  mostly  on  the  electrolyte 
matrix  side  of  the  electrode,  which  is  also  an- 
other inefficient  use  of  the  electrodes. 

Considering  the  previous  two  extreme  cases, 
there  ought  to  be  a  certain  intermediate  elec- 
trode composition  in  which  there  is  enough  cata- 
lyst to  allow  a  continuous  electrolyte  phase  and 
continuous  paths  for  electron  conduction,  and 
enough  Teflon  so  that  there  will  be  adequate  con- 
tinuous paths  for  the  transfer  of  gaseous  reac- 
tants.   Industrially,  electrodes  with  80  wt.  % 
Pt-black  have  been  considered  near  optimal  and 
are  most  widely  used,  but  no  entirely  adequate 
modeling  of  these  electrodes  has  yet  been 
achieved.   There  exists  a  need  to  identify  all 
transport  processes  occurring  in  operating  elec- 
trodes and  to  devise  a  mathematical  model  of  an 
electrode  which  will  make  it  possible  to  correlate 
known  electrode  behavior  with  the  predictions  of 
the  model,  and  to  design  improved  electrodes. 


CHAPTER  III 

ELECTRODE  REACTIONS,  IONIC  HYDRATION, 
AND  WATER  TRANSPORT 


3. 1   Introduction 

Alkaline  fuel  cells,  as  pointed  out  by  Kordesch  [6], 
are  technologically  the  most  advanced  fuel  cell  type  today. 
Despite  their  high  power  output  and  increasingly  lower  cost, 
many  problems  remain  to  be  solved  before  they  can  have  ex- 
tensive commercial  applications.   The  concentration  gradients 
of  electrolyte  in  the  matrix  of  an  operating  fuel  cell  have 
been  measured  by  Miller  and  Fornasar  [7],  and  by  Lundquist 
and  Vogel  [8].   However,  the  transport  of  free  water  within 
the  electrolyte  matrix  has  not  been  discussed  in  the  litera- 
ture.  In  this  chapter  we  will  look  into  the  problem  of 
water  transport  in  the  electrolyte  matrix  in  relation  to 
electrode  reactions  and  ionic  hydration  to  develop  ways  of 
estimating  the  concentration  gradients  of  free  water  and 
the  electrolyte  in  the  matrix. 

The  overall  reaction  product  in  the  alkaline  hydrogen- 
oxygen  fuel  cell  is  water.   Water  generated  by  this  reaction 
must  be  removed  from  the  cell  to  avoid  dilution  of  the  elec- 
trolyte, but  water  transport  within  the  fuel  cell  is  also 
very  important  in  another  sense:   it  is  both  a  reactant  at 
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the  cathode  and  a  solvator  of  the  hydroxyl  ions  formed 
there.   The  reactions  at  the  electrodes  can  be  written  as 
follows: 


Anode  reaction:   2  H^  +  4 [OH  -s  HO]  =  4  (l  +  s)  H^O  +  4e 
Cathode  reaction:   O-  +  2(l+2s)  H2O  +  4e  =  4[0H"''S  H^O] 
Overall  reaction:        2H2  +  0^  =  2H2O 


where  s  is  the  primary  hydration  number  of  the  hydroxyl 
ion. 

It  can  be  seen  that  water  must  be  transferred  to  the 
cathode  both  to  maintain  the  cathode  reaction  and  to  sol- 
vate the  hydroxyl  ions  produced  there.   Water  transport 
from  the  anode  to  the  cathode  involves  transport  of  sol- 
vent water  whereas  water  transport  from  the  cathode  to  the 
anode  involves  transport  of  hydrated  hydroxyl  ions. 

All  of  the  ions  in  the  electrolyte  are  hydrated  but 
our  principal  concern  is  with  the  hydroxyl  ions  as  these 
carry  most  of  the  intracell  current  and  they  participate 
in  the  electrode  reactions,  whereas  potassium  ions  do  not. 
As  pointed  out  by  Bockris  [9],  water  may  be  regarded  as 
being  bound  to  ions  in  two  ways:   (i)  strongly  bound  and 
(ii)  weakly  bound.   Strongly  bound  water  -  called  the  pri- 
mary hydration  sheath  -  is  held  so  strongly  to  ion  that  it 
moves  with  the  ion.   Weakly  bound  water  -  called  the 
secondary  hydration  sheath  -  is  attached  so  loosely  to  the 
ions  that  it  behaves  like  a  diffuse  cloud  around  the  ion 
and  is  in  continuous  interchange  with  water  in  the  surrounding 
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medium.   Only  the  primary  hydration  sheath  will  be  con- 
sidered here.   The  number  of  water  molecules  in  the  pri- 
mary hydration  sheath  of  an  ion  is  called  the  primary 
hydration  number  of  the  ion. 

3. 2   Mass  Transfer  in  the  Electrolyte  Matrix 

The  electrolyte  matrix  usually  consists  of  a  thin 
sheet  of  specially  processed  asbestos  with  the  voids  be- 
tween the  asbestos  fibers  being  filled  with  electrolyte, 
which,  in  this  case,  is  concentrated  potassium  hydroxide. 
All  transport  takes  place  in  the  electrolyte  solution 
which  has  three  constituents:   water,  potassium,  and  hydro- 
xyl  ions.   However,  these  ions  do  not  exist  as  bare  ions: 
rather,  they  exist  as  solvated  ions  and  the  water  in  the 
primary  hydration  sheath  around  the  ion  is  fairly  tightly 
bound  so  that  it  moves  with  the  ion. 

As  the  hydrated  hydroxyl  ion  moves  from  a  region  of 
high  KOH  concentration  to  one  of  low  KOH  concentration, 
the  hydration  number  increases  and  a  relationship  of  the 
following  form  can  be  written: 


OH  •s^H20  +  (s2-s-j^)  H2O  =  OH  •S2H2O 


Since  the  KOH  concentration  in  a  fuel  cell  is  usually  large 
and  since  according  to  Bockris  and  Saluja  [10]  the  hydration 
number  appears  to  change  slowly  with  KOH  concentration  in 
this  concentration  range,  it  can  be  assumed  that  s^  -  s-,  is 
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small  so  that  an  average  value  s  can  be  used  over  the 
entire  electrolyte  matrix. 

In  a  multicomponent  system  the  momentum  flux  depends 
only  upon  the  velocity  gradients  and  the  energy  flux 
depends  upon  the  temperature  gradient.   In  the  electrolyte 
matrix  considered,  although  there  may  be  some  local  varia- 
tions in  temperature,  the  temperature  differences  across 
the  matrix  are  small  compared  to  the  absolute  temperature. 
Therefore  constant  temperature  will  be  assumed  in  the 
electrolyte  matrix.   Also,  because  the  electrolyte  is  con- 
tained in  the  small  void  spaces  of  the  electrolyte  matrix, 
convection  will  be  assumed  to  be  absent,  i.e.,  the  bulk 
velocity  of  the  fluid  in  the  matrix  is  zero.   Hence,  it  is 
not  necessary  to  consider  the  equations  of  motion  and  energy. 

The  mass  or  molar  flux  of  a  species  depends  both  on 
mechanical  driving  forces  (ordinary,  pressure,  and  forced 
diffusion)  and  on  thermal  diffusion  effects.   The  multicom- 
ponent diffusion  equation  is  [11]: 
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Vy .  +  s . VT  -  —i  7P 
.   1     1       P 


=  RT  I 
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:  Vij  - 
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T     T 

D.    Dt 

1     1 


V2.nT 


(3.1-1) 


where   p.  -    electrochemical  potential  of  species  i.  Joule/mole; 

s.  =  partial  entropy  of  species  i,  Joule/mole-deg; 

2 
V .  ■    =   diffusion  coefficient  of  i  in  j,  cm  /sec; 

T 
D.  =  thermal  diffusion  coefficient  of  i,  gm/cm-sec; 

C   =  total  solution  concentration,  mole/1; 

V.  -    average  velocity  of  species  i,  cm/sec. 
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It  was  already  assumed  that  the  temperature  in  the 
electrolyte  matrix  is  constant  and  the  gradient  of  pres- 
sure is  usually  very  small.   Thus,  the  thermal  diffusion 
and  the  pressure  diffusion  terms  are  usually  of  only 
secondary  importance  and  may  be  neglected.   Equation  (3.3-1) 
then  takes  the  form 


C  .  C  . 

C.  Vy.  =  RT  I    ^y-  (V  -v^) 


(3.3-2) 


The  V . .    are  the  diffusion  coefficients  describing  the 
interactions  between  species  i  and  j.   For  potassium  hydro- 
xide solutions  we  have  two  independent  equations 
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C   Vy   =  RT 
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(3.3-4) 


where   C   =  concentration  of  total  free  water,  moles/1; 

C   =  concentration  of  potassium  ion,  moles/1; 

C_  =  concentration  of  hydroxyl  ion,  moles/1; 

C   =  total  solution  concentration,  moles/1; 

y   =  electrochemical  potential  of  the  potassium  ion. 
Joule/mole; 

y_  -   electrochemical  potential  of  the  hydroxyl  ion. 
Joule/mole. 


The  flux  equations  for  the  potassium  ion  and  the  hydrox- 
yl  ion  may  be  obtained  by  rearranging  Equations  (3.3-3)  and 


(3.3-4).   For  steady-state  conditions,  these  flux  equations 
can  be  written  as  follows  (see  Appendix  1  for  details) : 
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(3.3-6; 


where  V    is  the  diffusion  coefficient  of  the  electrolyte  and 
is  given  by 


V   = 


(3.3-7: 


and  V   is  the  diffusion  velocity  of  water.   The  transference 
o  -^ 

numbers  (with  respect  to  the  solvent  velocity)  are  defined 
as  follows: 


t.  =  1  -  t   =    ri  n 


(3.3-8) 


The  measured  diffusion  coefficient,  D,  is  related  to  V   by 


D  = 


VC       r-  d  £n  y^-. 


o  •- 


1  + 


d  £n  m 


(3.3-9) 


where  y+  is  the  mean  molal  activity  coefficient  and  m  is  the 
molality. 

It  can  be  shown  that  the  relation  between  the  gradient 
of  chemical  potential  and  the  gradient  of  concentration  is 
given  by  (see  Appendix  2) 


39 


o 


1  - 


d  £n  C  -I 
o 

d  Jin  C 


VC 


(3.3-10) 


Substitution  of  Equation  (3.3-10)  into  Equations  (3.3-5)  and 
(3,3-6)  gives 
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(3.3-12) 


For  water,  the  total  flux  consists  of  the  flux  due  to  free 
water,  that  due  to  water  of  hydration,  and  that  due  to  con- 
vection.  Since  potassium  ions  do  not  enter  into  the  elec- 
trode reactions,  the  net  flux  of  potassium  ions  is  zero  and 
the  flux  of  hydration  water  with  potassium  ions  is  zero. 
Thus,  we  have 


N„  ^  =  -DVC^   +  sN_  +  C  V 
-H_0        fw     -      o-o 


(3.3-13) 


where  C_   is  concentration  of  free  water, 
fw 

It  was  already  assumed  that  all  transport  is  essentially 
one-dimensional,  so  the  flux  equations  become 
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VC  +  — -  +   C_v 
z  F      -o 


(3.3-14) 


dC 


^H^O 


-D 


fw 


dx 


+  sN- 


(3.3-15) 
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In  order  to  solve  these  equations,  it  is  necessary  to  know 
the  transference  number  and  the  hydration  number  of  ions. 

3. 3   Transference  Numbers  in  KOH  Solutions 

The  transference  number  of  an  ion  is  the  fraction  of 
the  total  current  carried  by  that  ion.   Several  methods  for 
measuring  transference  numbers  are  known  and  several  em- 
pirical equations  for  transference  number  as  a  function  of 
concentration  have  been  developed.   However,  most  of  the 
aforementioned  studies  have  been  done  intensively  on  dilute 
solutions  and  studies  in  concentrated  electrolytic  solutions 
have  been  less  frequent.   Among  the  empirical  equations  is 
the  Onsager  equation,  according  to  which  the  equivalent  con- 
ductance of  an  ion.  A.,  is  given  by 


X.  =  X°  -   i-  A  +  BA° 
1     1     2        1 


/C 


(3.3-16) 


where  X.  -    ionic  conductance  of  species  i  at  infinite  dilu- 
tion, cm2/ohm; 

A,B  =  constants. 


Thus,  the  transference  number  of  the  potassium  ion  is  given 
by 


X^       X°  -  (|  A  +  BX°)/C 
'^    '  V^X  "  A°  -  (  A   +  BA°)/C 


(3.3-17) 


where  A  -    ^.    '^   ^_ 
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While  the  Onsager  theory  predicts  the  concentration  depend- 
ence of  the  transference  number,  it  is  applicable  only  at 
low  concentrations. 

Recently,  Merenkov  [12]  measured  the  transference  num- 
ber of  potassium  in  KOH  solutions  of  high  concentrations. 
It  was  found  that  the  transference  number  of  potassium  ion 
was  independent  of  temperature  from  18  to  65  °C  but  that  it 
decreased  with  increase  in  KOH  concentration.   The  transfer- 
ence number  of  potassium  ion  was  correlated  by  the  empirical 
equation 

t_^  =  0.260  -  0.047(/C  +  1)  (3.3-18) 

On  the  other  hand,  Knobel  [13]  determined  the  transfer- 
ence number  of  the  potassium  ion  in  KOH  solutions  by  the  EMF 
method;  he  found  that  the  transference  number  of  potassium 
was  0.2633  and  that  it  was  independent  of  concentration  over 
the  concentration  range  from  0.03  to  3.0  N.  The  predictions 
of  the  Onsager  theory  and  the  experimental  results  of  Knobel 
and  Merenkov  are  summarized  in  Table  3.1. 

In  fuel  cell  operations  concentrated  solutions  of  KOH 
are  used.   Therefore,  the  Merenkov' s  data,  which  are  the 
only  ones  available  for  high  concentrations,  will  be  used  in 
the  later  calculations. 


3 . 4   The  Primary  Hydration  Numbe r s 
in  Concentrated  KOH  Solutions 


The  hydration  number  of  an  ion  is  the  number  of  water 
molecules  attached  to  the  ion,  and  the  absolute  value  of  the 
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Table  3.1 
Summary  of  Transference  Number  Data  for  KOH  Solutions. 


Concentration  of 

KOH  Solutions  (N)   Onsager  Eq, 


Transference  Number  of  K   from 
Knobel      Merenkov 


0.03 
0.1 
0.3 
1.0 
3.0 
5.0 
7.0 
9.  0 
10.0 


0.261 
0.251 
0.232 
0.176 


0. 

2633 

0.205 

0. 

,2633 

0.198 

0, 

,2633 

0.187 

0. 

,2633 

0.166 

0. 

,2633 

0.132 

- 

0.108 

- 

0.088 

- 

0.072 

_ 

0.  064 

hydration  number  has  an  important  bearing  on  the  amount  of 
"free"  water  present  in  the  electrolytic  solutions.   For 
example,  in  30  wt.  %  KOH  solution,  if  the  hydration  number 
of  both  the  potassium  and  hydroxyl  ions  is  equal  to  one, 
then  the  concentration  of  "free"  water  is  36.4  moles/1,  while 
if  the  hydration  number  of  the  ions  is  equal  to  four,  there 
will  not  be  any  "free"  water  available.   The  concentration 
of  "free"  water  for  different  values  of  the  hydration  number 
at  various  KOH  concentrations  is  tabulated  in  Table  3.2 
(assuming  that  the  hydration  number  of  the  potassium  ion  is 
equal  to  that  of  the  hydroxyl  ion)  and  shown  in  Figure  3.1 
as  a  function  of  KOH  concentration. 
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Table  3.2 


Concentration  of  "Free"  Water 
as  a  Function  of  Hydration  Number. 


Concentration       Concentration  of  Free  Water,  C^  ,  moles/1 

rw 


of  KOH  Solution 

wt.  %    Moles/1   s  =  0 


10       1.91    53.58    49.76    45.93    42.11   38.29   34.46 
20       4.16    51. SI    43.49    35.16    26.84   18.52   10.20 


30  6. S3  45.  2"  35  .6"  ll. :5 
40  9.86  45.87  26,14  6.41 
50      13.44    41.52    14.65 


The  concept  of  Lonic  hydration  has  been  investigated 
by  many  workers  and  several  experimental  methods  have  been 
used  for  determining  the  hydration  number  of  ions.   However, 
the  different  experimental  methods  give  contradictory  results 
because  of  a  lack  of  a  really  precise  definition  of  the 
hydration  number.   Recently ,  Bockris  [10]  used  the  compressibility 
method  to  determine  the  sum  of  the  hydration  numbers  of  ions 
in  solutions  of  halides.   By  combining  compressibility 
measurements  of  the  ionic  vibration  potential,  individual 
hydration  numbers  for  certain  cations  and  anions  were  estimated, 
According  to  these  results,  the  hydration  number  of  potassium 
ion  in  dilute  solutions  is  3.5  while  that  for  hydroxyl  ion 
can  be  estimated  to  be  equal  to  4 . 3 .   If  these  hydration 
numbers  remain  constant  in  concentrated  solutions,  as  pointed 
out  earlier,  there  would  not  be  any  "free"  water  available  in 
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Figure  3.1:   Concentration  of  "Free" Water   for  Different 
Values  of  Hydration  Number. 
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30  wt.  %  KOH  solution.   It  is   obvious,  therefore,  that 
the  hydration  number  must  vary  with  the  concentration  of 
the  solution,  but  there  is  little  experimental  information 
concerning   the  concentration  dependence  of  the  hydration 
numbers  of  the  potassium  and  hydroxyl  ions  at  high  con- 
centrations. 

From  some  of  the  limited  information  available  in  Bockris' 
paper   [10],   the    hydration  number  of  the  potassium  and 
hydroxyl  ions  may  be  estimated  by  extrapolating  the  data 
available.   This  procedure  results,  for  30  wt.  %  KOH  solu- 
tion, in  an  estimated  hydration  number  of  potassium  ion  of 
2.0  and  2.0  for  the  hydroxyl  ion.   These  values  are,  however, 
very  tentative  owing  to  the  extrapolation  required.   Since 
they  are  the  best  estimates  available,  they  are  used  in  the 
next  section  to  estimate  concentration  gradient  in  the 
electrolyte  matrix  for  different  current  densities. 

3 . 5   Concentration  Gradients  in  the  Electrolyte  Matrix 

Having  estimated  the  hydration  numbers  and  the  trans- 
ference numbers  of  the  ions,  the  flux  equations  for  the 
potassium  ion,  the  hydroxyl  ion,  and  water  can  be  solved. 
The  concentration  gradients  of  the  "free"  water  and  that  of 
the  electrolyte  in  the  electrolyte  matrix  can  be  estimated 
as  follows: 

3.5.1   Concentration  Gradients  of  "Free"  Water 

For  the  "free"  water  concentration  gradients  in  the  x- 
direction,  we  have  from  Equation  (3.3-15) 
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dC 


From  the  electrode  reactions  at  both  anode  and  cathode  it 

can  be  seen  that  the  flux  of  "free"  water  from  the  anode  to 

the  cathode  is  one-half  of  the  hydroxyl  ion  flux  in  opposite 

direction.   Therefore,  the  following  can  be  written. 

dC 
^H20,x  =  -"^  -d^  ^  ^N-,x  =  -^-,x/2         (3.5-1) 

Since,  according  to  the  Faraday's  law,  N    =  -^  ,  the  "free" 

—  ,  X    Z  r 
water  concentration  gradient  can  be  expressed  in  terms  of  the 

hydration  number  of  ions,  the  diffusivity,  and  the  current 

density,  that  is 


dC^ 

fw 

dx 


=  ^l 


z  DF 


(3.5-2; 


Values  of  the  diffusion  coefficient  of  the  electrolyte  at 
various  concentrations,  taken  from  Bhatia's  results  [14], 
are  converted  into  the  molar  average  frame  (see  Appendix  3) . 
With  the  available  data  on  diffusivity  and  the  various  values 
of  the  hydration  number,  the  concentration  gradients  of  the 
"free"  water  at  different  KOH  concentrations  can  be  calculated 
as  a  function  of  current  density.   The  results  of  "free" 
water  concentration  gradients  are  summarized  in  Table  3.3. 

3.5.2   Concentration  Gradients  of  KOH  Solution 

For  the  concentration  gradient  of  the  KOH  solution,  we 
have  from  Equation  (3.3-14) 
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Table  3.3 
"Free"  Water  Concentration  Gradients, 


Hydration  No.  of 
Hydroxyl  ion,  s 


KOH  Concentration,  C,  moles/1 
1.91      4.16      6.80      9.86      13.44 


0 
1 
2 
3 
4 
5 


0. 085i 
0.236i 
0.365i 
0.470i 
0.552i 
0.610i 


0. 078i 
0.199i 
0.274i 
0. 303i 
0.285i 
0.220i 


0.075i 
0.171i 
0.221i 
0.202i 


0.074i 
0.143i 
O.lOBi 


0.071i 
O.llOi 


N     =  -vD 
-,x 


d  £n  C 


1  - 


fw 


dC 


it 


+ 


+  C  V 


d  Jin  C   J  dx    z_F     -  o,x 


(3.5-3) 


N 


Since  N     =  — -    ,  1  -  t   =  t, ,  and  N„  ^ 

-,x    z_F         -     +        H„0,x 

Equation  (3.5-3)  yields 


-^   =  C.  V 
2      f  w  o ,  X 


dC 
dx 


The  term 


t .  + 


z  F    +     2C 


fw 


/v  D 


d  £n  C 


1  - 


fw 


d  £n  C 


(3.5-4) 


d  e.n  C 


1  - 


fw 


d  en  C 


can  be  shown  to  be  equal  to 


where  V   is  the  partial  molal  volume  of  the  solvent 

o         '^ 

water  and  is  given  by 


C.   V 
fw  o 


M 


V   = 


o 


P  -  C 


dp 
dC 


(3.5-5) 


where  M   =  molecular  weight  of  water, 
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Therefore,  Equation  (3.5-4)  becomes 


dC 

C.   M 
fw   o 

i^^  '  -.J 

dx 

v_z_DF 

(-^^) 

(3.5-6] 


As  for  the  density-concentration  relation  of  KOH  solution, 
an  equation  that  relates  the  density  with  the  concentration 
has  been  developed  from  the  density-concentration  data [15] by 
using  the  Lagrange  interpolation  formula. 

p(C)  =  1.00086  +  C(0. 04930  -  C(0.0143  -  C(0. 000059 

-  0.0000014C) ) )  (3.5-7) 

Differentiation  of  Equation  (3.5-7)  gives 

^  =  0.04930  -  C(0. 00286  -  C(0. 000177  -  0.0000056C)) 

(3.5-8) 
The  concentration  gradients  of  KOH  solution  calculated  from 
Equation  (3.5-6)  are  summarized  in  Table  3.4. 

3 . 6   Discussion 


From  Table  3.4  it  can  be  seen  that,  for  the  estimated 

values  of  the  transference  and  hydration  numbers,  a  current 

2 
density  of  one  ampere  per  cm  ,  and  a  mean  KOH  concentration 

of  30  wt.  %  (=  6.80  N) ,  the  KOH  concentration  gradient  in 

the  electrolyte  matrix  is  0.009  gmole/1-cm.   Since  the 

electrolyte  matrix  is  normally  less  than  one  millimeter  in 

thickness,  it  seems  clear  that  the  KOH  concentration  dif- 
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Table  3.4 
KOH  Concentration  Gradients  in  Electrolyte  Matrix, 


KOH  Concentration,  C,  moles/1 
Hydration  No.  of 
Hydroxyl  Ion,  s     c  =  1.91      4.16        6.80        9.86 


1  0. 02431  0. 01721     0.01331     0. 01091 

2  0. 02101  0.01231     0.00901     0.00401 

3  0.01791  0.00821     0. 00481 

4  0.01511  0.00491 

5  0.01251  0.00241 

ference  between  the  anode  and  the  cathode  even  at  the  high 

2 
current  density  of  one  ampere  per  cm   appears  to  be  almost 

negligible. 


CHAPTER  IV 

MATHEMATICAL  MODEL  OF  TEFLON-CATALYST  LAYER  AND 
CALCULATION  OF  CURRENT  DENSITY  DISTRIBUTION 


4 . 1   Introduction 

Investigation  of  the  electrode  structure  by  scanning 
electron  microscopy  (SEM)  and  energy  dispersive  X-ray 
analysis  (EDXA)  has  yielded  much  information  on  how  the 
platinum  black  catalyst  and  the  Teflon  particles  are  mixed 
to  form  the  catalyst-Teflon  layer  of  the  electrode.   On 
the  basis  of  this  information  we  can  understand  how  the 
electrode  is  wetted,  how  gaseous  reactants  can  be  trans- 
ported, and  what  the  means  of  ionic,  electronic,  and  water 
transport  are.   Mathematical  modeling  of  these  catalyst- 
Teflon  layers  has  been  done  to  relate  these  findings  to 
predicted  performance  of  the  electrodes. 

A  scanning  electron  micrograph  of  the  surface  of  a 
layer  of  Pt:Teflon  in  the  weight  ratio  80:20  is  shown  in 
Figure  4.1  (Magnification  =  72,000).   From  studies  of  this 
and  similar  SEM  the  following  conclusions  can  be  drawn: 

1)   The  Teflon  aifqi'egates  are  prolate  spheroids  of 
sizes  2000  to  3000  A  in  the  major  axis  and  1000 
to  1500  A  in  the  minor  axis.   However,  the  sub- 
structure of  these  spheroids  shows  that  they 
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consist  of  aggregates  of  spherical  Teflon  mole- 
cules with  diameters  of  200  to  300  A. 


Figure  4.1:   Scanning  Electron  Micrograph  of  the  Surface 
of  the  80:20  Pt-black-Tef Ion  Layer, 
Magnification:   72,000  . 


2)  Owing  to  thin  extremely  small  size  the  catalyst 
crytallites  cannot  be  resolved  in  the  SEM,  but 

it  is  clear  that  they  have  an  equivalent  diameter 
of  less  than  100  A.   From  measurement  of  its  sur- 
face area  by  the  BET  gas  adsorption  method  and 
the  estimation  of  its  packing  condition  as  shown 
in  Chapter  II,  it  seems  very  likely  that  the  dia- 
meter of  a  single  catalyst  crystallite  is  in  the 
range  of  30  to  40  A. 

3)  The  Teflon  aggregates  are  coated  with  thin  layers 
of  catalyst  crystallites  and  the  distribution  seems 
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to  be  relatively  uniform.   As  the  ratio  of 
platinum  black  to  Teflon  increases,  the  chances 
of  observing  bare  Teflon  decrease,  and  the  great 
majority  of  Teflon  aggregates  are  coated  by 
platinum  black  in  the  electrode  of  80:20  platinum 
black-Teflon  mixture 

From  these  observations  it  appears  that  the  following 
processes  occur  during  the  formation  of  an  electrode: 

1)  When  a  well-dispersed  suspension  of  platinum  black 
is  mixed  with  Teflon  30  dispersion  (which  is  nega- 
tively charged  hydrophobic  colloid) ,  the  Teflon 
aggregates  are  coated  with  platinum  black. 

2)  During  the  filtration  of  this  mixture  the  catalyst- 
coated  aggregates  settle  on  the  filter  medium  and 
form  a  layer  of  platinum  black-coated  Teflon  par- 
ticles. 

3)  When  the  filter  cake  is  pressed,  these  particles 
are  compacted  together  very  tightly. 

4)  During  the  sintering  process,  the  Teflon  aggregates 
expand,  soften,  and  contact  one  another  to  form  a 
bonded  layer  with  more  or  less  continuous  structure, 

On  the  basis  of  these  observations  and  conjectures  the 
electrode  can  be  pictured  as  consisting  of  gas  channels  of 
varying  cross  section  through  the  intra-aggregate  voids  of 
the  Teflon  aggregates  which  are  themselves  coated  with  a 
catalyst  layer  on  their  external  surface.   This  catalyst 
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layer    is   wetted    and   penetrated    by    the    electrolyte.       Thus, 
gas   penetrates    the    catalyst-Teflon    layer   by   means    of    the 
intra-aggregate    voids    and   electrolyte    penetrates    the    extra- 
aggregate    voids   where    the    catalyst   crystallites    lie.       Gas, 
electrolyte,    and   catalyst   meet    at    the    surface    of    the    Teflon 
aggregates    as    shown    schematically    in    cross    section    in 
Figure    4.2. 


CATALYST 
CRYSTALLITES 


ELECTROLYTE 
MATRIX 


TEFLON 
AGGREGATE 

TEFLGN 
rOLECULE 

INTRA- 
AGGREGATE 
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VOIDS 

CURRENT 
COLLECTOR 
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Figure  4.2:   Schematic  Diagram  of  Platinum  Black  Catalyst 
and  Teflon  Particles  in  Electrodes. 


A  cylindrical  multicompartment  model  as  shown  in 
Figure  4.3  is  a  much  simplified  but  fairly  reasonable  approx- 
imation of  the  structure  of  the  catalyst-Teflon  layer.   A 
compartment  in  the  model  represents  a  Teflon  aggregate  which 
consists  of  several  Teflon  molecules  and  interconnected  intra- 
aggregate  voids.   Windows  represent  narrow  gas-phase  connec- 
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tions  between  Teflon  aggregates  which  may  limit  the  transport 
of  gaseous  reactants. 


TEFLON  WITH 
INTRA-AGGREGATE  VOIDS 


i 
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Figure  4.3:   Schematic  Model  of  the  Structure  of  the 
Catalyst-Teflon  Layer. 


Although  the  model  described  above  represents  the  struc- 
ture of  electrode  layer  very  well,  there  arises  a  problem  in 
mathematical  modeling  due  to  the  existence  of  the  windows 
which  limit  the  transport  of  gaseous  reactants.   For  the 
sake  of  simplicity  in  mathematical  modeling  and  due  to  the 
limitation  in  time,  the  model  shown  in  Figure  4.4,  which  is 
the  same  as  the  multicompartment  model  except  having  no 
windows,  is  adopted  for  this  work. 

For  an  electrode  modeling  only  the  oxygen  electrode  was 
considered  because  the  limiting  current  density  of  a  hydrogen- 
oxygen  fuel  cell  is   usually  determined  by  the  oxygen  electrode 
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Figure  4.4:   Schematic  Model  of  the  Structure  of  the 
Catalyst-Teflon  Layer  with  No  Windows. 


owing  to  the  slowness  of  the  reduction  of  oxygen  on  the 
catalyst  surface.   The  following  assumptions  were  made  in 
the  development  of  the  mathematical  model  to  simplify  the 
computations : 


1)  Steady  state,  i.e.,  constant  current  density. 

2)  Isothermal. 

3)  Diffusion  coefficients  and  other  physical  pro- 
perties are  independent  of  electrolyte  concentra- 
tion. 

4)  Potential  and  current  are  assumed  to  vary  in  the 
x-direction  only. 
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5)  The  concentration  of  KOH  is  assumed  to  vary  in 
the  x-direction  only. 

6)  The  diffusion  of  oxygen  through  gaseous  space  is 
not  restricted. 

7)  Both  cations  and  anions  have  same  primary  hydration 
number.   In  this  work  it  is  taken  to  be  2. 

4 . 2   Electrode  Reaction  and  Reaction  Rate 

The  electrochemical  reaction  in  the  oxygen  electrode 
may  be  written  as 

O2  +  10  H2O  +  4  e~  ^  4[OH~'2H20] 

For  convenience  we  denote  OH  "2^120,  K  •2HpO,  O2  /  and  H^O 
as  species  1,  2,  3,  and  4,  respectively,  in  the  equations 
which  follow. 

Although  several  workers  [16,17]  have  studied  the 
mechanism  of  the  reduction  of  oxygen  in  alkaline  solutions 
on  platinum  catalysts,  it  is  still  inadequately  understood. 
However,  the  work  of  Gnanamuthu  and  Petrocelli  [18]  suggested 
that  a  suitable  rate  equation  is 


di 

c 

dx" 


=  ai 


C^(x)-^  (C^ix)-^ 


,0 


exp(-F(4,^-(},^)/RT) 


C^(x) 


c. 


exp(F((})g-(})5,)/RT) 


(4.2-1) 


57 

2 
where   i   =  current  density  in  the  catalyst  phase,  amp/cm  ; 

2    3 
a   =  surface  area  per  unit  volume  of  electrode,  cm  /cm  ; 

2 

1   =  exchange  current  density,  amp/cm  ; 

C.  =  concentration  of  species  i,  gmole/1; 

C.  =  concentration  of  species  i  at  x  =  0,  gmole/1; 

(J)   =  potential  in  the  catalyst  phase,  volt; 

()).  =  potential  in  the  electrolyte  phase,  volt. 

Equation  4.2-1  must  be  solved  with  the  following  boundary 
conditions: 


At  X  =  0,      i   =0 

s 

At  X  =  L,      i   =  i 

s     app 


where  i     =  apparent  current  density  being  drawn  from  the 
app    ^'^  J  ^ 

2 

electrode,  amp/cm  . 


4 . 3   Ionic  Fluxes  in  the  Electrolyte  Phase 

At  steady  state  the  mass  fluxes  in  the  x-direction  for 
a  binary  electrolyte  system  can  be  written  as  [19] 


P*v,C  dp     i„t* 

N,    = ^  :^  +   -\t^  +  C,v*  (4.3-1) 

l,x       vRTdx     ZtF      11, X 


V*VyC    du     ij,t* 

N_    =  0  = -^  -^  +   — -^  +  C„v*        (4.3-2) 

2,x  vRTdx     z„F      22, X 


where  N..    and  N     are  fluxes  of  hydroxyl  and  potassium  ions 
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Since  C^  =  C^^  +  C2  +  C^  =  2C  +  C^  , 

44, x     4,x      l,x^' 
Equation  (4.3-2)  can  be  rearranged  to  give 


V*v^C    du 
2     e 

vRT   dx 


t* 


Cr 


(4.3-3) 


Since  ^a^  =  vRT2,na,  z^  =  1,  and  D* 


D 


V 


Appendix  3),  Equation  (4.3-3)  becomes 


(C+C4)V4 


^.    dJlna  , 
=  ^   dl^  (^^^ 


2 dg.nC 

dx 


(C+C^)V^ 


t*    -    ^ 
2    C^ 


(4.3-4) 


C+C4t2 
Substitution  of  t*  =  — into  Equation  (4.3-4)  gives 


dC  ^  ^R^2^4^4 
dx      ^^v^ 


(C  +  C^) 


(4.3-5) 


which  will  be  solved  by  using  one  of  the  boundary  conditions, 

namely,  at  x  =  0 ,  C  =  C  . 
^  o 

4.4   Potential  Gradients  in  the  Electrode  Layer 


The  potential  drop  in  the  catalyst  phase  is  governed  by 

Ohm' s  law 

d(t) 

(4.4-1) 


1   =  -o  -2 — 
s       dx 
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where  a  is  the  conductivity  of  catalyst. 

The  current  carried  by  the  ions  is  related  to  the  poten- 
tial (})   by  the  following  equation  [20] 


d4), 


i   =  K  — - 
I  dx 


F 


-4^ 


nv. 


^2^2 


nc 


dx 


(4.4-2; 


where  k  is  the  conductivity  of  the  electrolyte  and  s.  is 
the  stoichiometric  coefficient  of  species  i  in  the  equation 


^1       ^2 
s,M,   +  s^M^   +  s  .M.  -*■   ne 


Substituting  yx      =   vRTJina  into  Equation  (4.4-2)  for  the  elec- 


trode reaction  equation  gives 


d<J), 


1         -1 
app      s 


dx 


VRT 
nF 


t      +    2^^ 
2  C, 


dg.na   dC 
dJlnC    dx 


(4.4-3) 


Equations  (4.4-1)  and  (4.4-3)  can  be  solved  with  two  B.C.'s, 

i.e.,atx  =  0,  (|)   =0,  and  at  x  =  L,  i   =  i 

^s     '  '   s     app 


4 . 5   Oxygen  Transport  in  the  Electrolyte  Phase 


The  concentration  of  oxygen  in  the  electrolyte  phase 
will  vary  both  in  the  r-  and  x-directions .   In  work  by  Tham 
[21],  (j)  ,  (J)  ,  C,  ,  Q.^'    ^^^   C.  were  assumed  to  vary  in  the  x- 
direction  only,  and  these  led  to  predicted  current  densities 
which  agreed  reasonably  well  with  experimental  values,  i.e., 
radial  concentration  gradients  appeared  to  be  relatively 
small.   To  evaluate  the  variation  of  Q^{r,y.)    for  a  given  x, 
an  average  value  of  C^ (x)  was  evaluated  at  that  x. 
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The  flux  equation  for  oxygen  is 


d  C^    D_  dC^    -,      di 

D   ^  +  _3  _^  _  i s  ^  Q 

3  ,2     r  dr     nF  dx 


(4.5-1) 


di 
where  -3—^  is  given  by  Equation  (4.2-1). 

After  several  mathematical  operations  (see  Appendix  4),  the 
following  equation  was  obtained  for  C-,(x). 


o 


C^{x)         Cj^(x)C^ 


C°C^(x) 


exp 


r2F(*3-ct,)- 


RT 


2r 


+ 


q (DEN) 


C^(x)C° 

1  -  — exp 

C°C4(x) 


r2F(t3-*,)i 


RT 


x  (I^(q(r^+6))K^(qrQ)  -  K^ (q (r^+6 ) ) I^ (qr^) ) 


(2r  6  +  6' 
o 


(4.5-2) 


ai 


where  q   = 


C,(x)      r   F((^^-(t>Jn 


o     c° 
nFCp^    4 


exp 


s  ^l' 


RT 


I  ,K   =  Bessel  functions  of  zero  order; 
o   o 

I,,K,  =  Bessel  functions  of  first  order; 

DEN  =  K, (q(r^+6) ) I^(qr^)  +  K  (qr^) I, (q (r  +5 ) ) 


Substitution  of  the  average  value  of  C^  (x)  for  C^ (x)  into 
Equation  (4.2-1)  gives 


di  ,C^(x)  C . (x) 

S  .     I     3                4        „ 

-5 =  ai   -7TT5 exp 

dx  o   „o     C"; 


F(4>3-*,) 


RT 


C^(x) 


exp 


rF(*3-*,; 


RT 
(4.5-3) 


61 


To  solve  for  C  (x)  ,  C^  (x)  ,  i  (x)  ,  (j)  (x)  ,  and  (})  (x)  ,  Equations 
(4.3-5),  (4.4-1),  (4.4-3),  (4.5-2),  and  (4.5-3)  can  be  com- 
bined with  the  accompanying  boundary  conditions. 

4 . 6   Computation  of  the  Current  Generation 

Equations  (4.3-5),  (4.4-1),  (4.4-3),  (4.5-2)  and  (4.5-3) 
were  integrated  numerically  using  the  quadratic  Runge-Kutta 
method;  the  computer  program  for  the  computation  of  the  cur- 
rent density  distribution  is  given  in  Appendix  5.   Numerical 
values  of  the  properties  of  the  electrolyte  such  as  electro- 
lytic conductivity,  partial  molal  volume,  and  the  activity 
of  the  electrolyte,  and  other  necessary  terms  were  expressed 
in  the  form  of  polynomial  equations  as  functions  of  the  KOH 
concentration.   The  specific  conductance  of  KOH  was  taken  from 
the  data  by  Klochko  and  Godneva  [22]  and  was  expressed  in  a 
polynomial  form  by  using  the  Lagrange  interpolation  formula 
as  follows: 

k:  =  -0.59876  +  C(0. 608757  -  C(0. 0684246  -  0 .  00229369C)  ) 

(4.6-1) 

In  the  calculation  of  current  generation  the  specific  surface 

2 
area  of  Pt-black  catalyst  was  taken  to  be  45  m  /g  (from  Chap- 

—  8 
ter  II)  and  an  exchange  current  density  for  oxygen  of  5.0x10 

2 
amp/cm   was  used  [23]. 
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4 . 7   Results  and  Discussion 

The  calculated  fractional  current  generation  is  given 
in  Table  4,1  and  the  current  generation  profiles  in  the  elec- 
trode layer  are  shown  in  Figures  4.5  and  4.6  as  functions  of 
catalyst  layer  thickness  and  current  density.   These  results 
agree  reasonably  well  with  those  obtained  from  actual  operating 
fuel  cells  [24,25] . 

From  Figures  4.5  and  4.6  one  can  see  that  at  low  current 
densities  the  distribution  of  the  current  generation  is  quite 
uniform,  while  at  high  current  densities  the  majority  of  the 
current  is  generated  in  that  portion  of  the  catalyst-Teflon 
layer  near  the  electrolyte  matrix. 

Cell  potential  and  power  density  are  tabulated  in  Table 
4.2,  and  these  are  plotted  as  a  function  of  current  density 

in  Figure  4.7.   From  Figure  4.7  we  can  see  that  the  limiting 

2 
current  density  is  approximately  1.3  amp/cm  ,  and  the  maximum 

2 
power  density  (0.70  watts/cm  )  is  obtained  at  a  current  density 

2 
of  about  1.00  amp/cm  . 

It  should  be  noted,  however,  that  in  this  work  it  was 
assumed  that  there  were  no  gas  transport  limitations;  moreover, 
it  was  also  assumed  that  the  oxygen  contained  no  impurities 
which  could  build  up  in  the  gas  spaces  and  introduce  gaseous 
concentration  gradients.   Thus,  the  model  for  which  calcula- 
tions are  presented  probably  overestimates  to  some  degree  the 
proportion  of  the  current  generated  in  the  catalyst  layer  near 
the  electrolyte  matrix.   Nevertheless,  it  seems  reasonable 
that  at  higher  current  densities  most  of  the  current  will  be 
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Distance  in  the  Electrode  Layer,  cm 


Figure  4.5:   Current  Generation  Distribution  in  the  Elec- 
trode Layer. 
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Figure  4.6: 


Cumulative  Fraction  of  Current  Generation  in 
the  Electrode  Layer. 
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Table  4.2 
Cell  Potential  and  Power  Density  vs.  Current  Density. 


Current  Density,        Cell  Potential,        Povs'er, 
amp/cm2  volt  watt/cm^ 


0.010  1.085  0.011 

0.050  1.058  0.053 

0.100  1.039  0.104 

0.200  1.009  0.202 

0.300  0.979  0.294 

0.400  0.948  0.379 

0.500  0.914  0.457 

0.600  0.878  0.527 

0.700  0.841  0.589 

0.800  0.802  0.641 

0.900  0.762  0.686 

1.000  0.704  0.704 

1.100  0.631  0.694 

1.200  0.489  0.587 
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Figure  4.7:   Cell  Potential  and  Power  Density  vs.  Current 
Density. 
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generated  near  the  electrolyte  matrix  because  of  the 
enormous  difference  between  liquid-phase  and  gas-phase 
transport  coefficients,  and  because  of  the  low  solubility 
of  oxygen  in  concentrated  KOH. 


CHAPTER  V 
SUMMARY  AND  CONCLUSION 

After  a  brief  description  of  fuel  cells  and  their  per- 
formance as  these  are  related  to  the  thermodynamics  and 
kinetics  of  the  electrode  reactions,  the  transport  processes 
occurring  in  an  alkaline  hydrogen-oxygen  fuel  cell  are  dis- 
cussed along  with  the  problems  associated  with  the  design  of 
electrodes  having  improved  performance.   However,  the  prin- 
cipal studies  completed  by  the  author  and  reported  in  this 
dissertation  are  concerned  with:   1)  determination  of  the 
surface  areas  of  electrode  components  and  evaluation  of  the 
effect  of  processing  temperatures  on  surface  areas,  2)  estima- 
tion of  particle  sizes  and  geometry  of  the  components  in  the 
electrodes,  3)  investigation  of  the  details  of  the  micro- 
structure  of  fuel  cell  electrodes  by  means  of  scanning  elec- 
tron microscopy  and  energy  dispersive  X-ray  analysis  resulting 
in  the  development  of  a  new  physical  picture  of  electrode 
structure,  4)  evaluation  of  the  influence  of  ion  hydration 
on  concentration  gradients  and  mass  transport,  5)  development 
of  a  physical  and  mathematical  model  for  analysis,  6)  develop- 
ment of  a  computer  program  for  use  in  predicting  electrode 
behavior,  and  7)  computation  and  analysis  of  current  genera- 
tion distribution  in  the  electrode  layer. 
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Surface  areas  of  platinum  catalyst,  Teflon  30  particles, 
and  their  mixtures  in  different  ratios  were  measured  by  the 
BET  gas  adsorption  method;  these  have  been  compared  with 
geonetrical  calculations  to  provide  information  as  to  their 
sizes  and  the  microstructure  of  the  catalyst-Teflon  layer, 
particularly  in  relation  to  the  results  obtained  by  scanning 
electron  microscopy.   The  dependence  of  surface  area  of 
platinum  black  on  the  temperature  of  pretreatment  was  also 
measured  and  related  to  electrode  processing  procedures.   It 
was  also  found  that  the  surface  area  of  the  platinum  catalyst 
depends  strongly  on  the  pretreatment  temperature  range  of 
310  to  350  ^C.   This  finding  was  very  important  as  the  nominal 
sintering  temperature  during  electrode  fabrication  is  about 
310 °C.   For  example,  it  was  shown  that  present  electrode 
sintering  temperatures  are  very  close  to  the  temperature  at 
which  the  catalyst  itself  is  sintered,  producing  a  large  re- 
duction in  surface  area.   The  need  for  very  good  temperature 
control  during  electrode  sintering  was  demonstrated. 

Scanning  electron  microscopy  and  energy  dispersive  X-ray 
analysis  techniques  were  used  in  investigating  the  detailed 
structure  of  pure  Teflon  and  catalyst  layer  of  different 
platinum  black-Teflon  ratios.   This  provided  information  on 
the  size,  shape,  and  microstructure  of  the  catalyst-Teflon 
layer  and  on  the  arrangement  and  distribution  of  catalyst  and 
Teflon  aggregates.   These  studies  showed,  for  example,  that 
Teflon  30  particles  are  aggregates  of  smaller  particles  which 
appear  to  be  Teflon  molecules.   Separate  calculations,  based 
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on  the  mean  molecular  weight  of  the  Teflon  and  its  density, 
conform  the  conclusions  drawn  from  electron  microscopy.   The 
aggregates  have  a  prolate  spheroidal  shape  with  minor  and 
major  diameters  of  1000  to  1500  A  (100  to  150  nm)  and  2000 
to  3000  A  (200  to  300  nm) ,  respectively,  while  the  molecules 
appear  to  be  spheres  with  diameters  ranging  from  200  to  300 
A  (20  to  30  nm) .   In  a  commonly  used  electrode  composition 
of  80:20  platinum  black  and  Teflon,  the  Teflon  aggregates 
are  seen  to  be  coated  with  one  or  more  layers  of  catalyst 
crystallites.   This  is  the  exact  opposite  of  the  previously 
held  view  of  the  structure  of  the  catalyst-Teflon  layer, 
namely,  that  this  layer  consists  of  roughly  spherical  ag- 
glomerates of  catalyst  (diameter  about  one  micrometer) 
surrounded  by  prolate  spheroidal  Teflon  particles.   The  sub- 
structure of  the  Teflon  particles  has  not  been  previously 
reported. 

The  size  of  individual  catalyst  crystallites  could  not 
be  determined  clearly  by  scanning  electron  microscopy  owing 
to  the  resolution  limitations  of  the  technique,  but  these 
studies  did  enable  us  to  determine  that  the  crystallite  size 
was  less  than  100  A  (10  nm) .   However,  this  information, 
combined  with  that  derived  from  the  surface  area  studies, 
made  it  possible  to  demonstrate  that  the  catalyst  crystallite 
size  is  almost  certainly  in  the  range  of  30  to  50  A  ( 3  to  5 
nm)  . 

Ion  hydration,  which  has  been  ignored  in  previously 
published  work,  has  been  studied.   Since  water  is  a  reactant 
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at  the  oxygen  electrode  and  a  product  at  the  hydrogen  elec- 
trode and  since  its  concentration  influences  the  concentration 
of  electrolyte  components,  it  is  reasonable  to  expect  that 
ion  hydration  will  affect  electrode  behavior.   Indeed,  ion 
hydration  was  found  to  have  a  significant  influence  on  the 
concentration  of  "free"  water.   In  this  work  a  primary  hydra- 
tion number  of  two  for  both  potassium  and  hydroxyl  ions  seems 
reasonable  for  the  KOH  concentration  used  in  fuel  cells,  and 
this  value  was  used  in  the  calculation  of  the  mathematical 
modeling  equations.   Unfortunately,  limitations  of  time  made 
it  impossible  to  examine  in  detail  the  influence  of  ionic 
hydration  in  the  mathematical  modeling;  that  will  be  a  prin- 
cipal objective  of  any  further  research  in  this  area. 

Based  on  the  information  obtained  earlier  from  the  studies 
on  the  size  and  distribution  of  electrode  components  and  on 
the  structure  of  electrodes,  a  schematic  model  for  mathematical 
analysis  was  devised  to  utilize  this  information  in  calculating 
concentration  and  current  generation  profiles.   A  set  of 
modeling  equations  was  developed  and  computations  done.   The 
results  were  analyzed  to  provide  information  on  current 
generation  distribution  in  the  electrode  layer  for  different 
current  densities. 

The  distribution  of  current  generation  and  the  magnitude 
of  the  limiting  current  density  given  by  mathematical  modeling 
of  the  oxygen  electrode  show  that  at  higher  current  densities 
most  of  the  current  is  generated  in  the  region  near  the  elec- 
trolyte matrix.   It  is  also  clear  that  there  is  a  need  to 
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determine  more  adequately  the  effective  mass  transport 
properties  of  the  electrode  layer  and  to  develop  and  analyze 
more  realistic  models.   Considering  the  fact  that  the  intra- 
aggregate  voids  of  Teflon  particles  have  very  narrow,  tortuous, 
and  interconnected  gaseous  paths  and  that  the  electrolyte 
wets  the  catalyst  layer  coated  on  the  surface  of  the  Teflon 
aggregates,  it  seems  very  likely  that  there  could  be  some 
hindrance  to  mass  transport  of  oxygen  in  the  catalyst-Teflon 
layer.   This  point  should  be  investigated  in  future  work, 
as  should  the  influence  of  impurities  in  the  oxygen  system. 


APPENDIX  1 

DERIVATION  OF  FLUX  EQUATIONS 
(EQUATIONS  (3.3-5)  AND  (3.3-6)) 


From  Faraday's  Law,  we  have 


i  =  F  y  z .N. 


X  1 


=  F  y  z .C. V. 

V    111 


(A. 1-1) 


Since  the  solution,  as  a  whole,  is  electrically  neutral, 

we  have 

y  z.C.  =  0  (A. 1-2) 

.11 

1 

Thus  Equation  (A. 1-1)  can  be  written  as 


=  F  y  z.C.v.  -  F  y   z.C.v 

V  1  1  1  h         1  1   o 
1  1 

-  F  y  z.C.   (v.  -  V  ) 

V  1  1     1  o 


(A. 1-3) 


Therefore,  for  a  1-1  electrolyte  such  as  KOH,  we  have 


i  =  Fz^C^  (v^  -  V  )  +  Fz  C   (v   -  v^; 
+  ++     o       ___     o 


From  Equations  (3.3-3)  and  (3.3-4),  we  have 


(A. 1-4) 


C^Vm^ 


=  RT 


+  o   - 


c_c^ 


(V 


'T  +- 


-  ^+) 


(3.3-3) 
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C  Vy   =  RT 


C_C  C_C^ 


)       (3.3-4; 


T  o- 


'T  +- 


Adding  the  last  two  equations,  we  get 


RTC  C  RTC_C 

C^Vy^  4-  C_Vy_  =  -^  (v^  -  v^)  +  ^-^    (v^  -  v_ )  (A.  1-5) 

T  o+  T  o- 


The  electrochemical  potential  of  the  electrolyte,  y  ,  is 
given  by 


Ue  =  ^+^+  +  ^-^. 


(A. 1-6) 


The  stoichiometric  concentration  of  the  electrolyte  can  be 
written  as 

C=-=^  (A.1-7) 


Hence 


C+y+  +  C_u_    =    Civ^M^    +    v_y_) 


=  Cy. 


Therefore,  we  get 


RTC  C 
CVy^  =  7^-,T — -    (v„  -  v^)  + 


RTC  C 


o 


e  -  C„0     -o    V'  '  C„P     ^^o  -  ^-^     ^^-^-^^ 


'T  0+ 


"T    o- 


Eliminating  v   from  Equations  (A. 1-4)  and  (A. 1-8),  we  get 


RTC  C  i 

-C^^   +  Fz^C^CVy^ 
T  0+ 


RTC  C  RTC_C 

^-p-°  Fz_C_  (v_  -  v^)  +  --^  F-+C^(v^  -  v_) 


'T  0  + 


'T  O- 
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RTC  C  C  F 
o  +  - 


r    z 


V    ^  V 

o+  o- 


(v^  -  vj 


RTC  C , C  F 
o  +  - 


zV   ^    -    z    V 
+   o+ -  o- 

o+  o- 


(v^  -  v_) 


(A. 1-9) 


Therefore 

iz  V  C^CVy  zV      V 

^     .         _         s     ^    -  o- T   "^0  +  0+  o- 

^_^v     v_;    ^  Y{zV   ^   -    z    V  RTC  (z  0  ^  -  z  P   ) 
-    +  o+     -  o-       O   +  0+     -  o- 

(A.1-10) 
Substituting  Equations  (3.3-7)  and  (3.3-8)  into  Equations 
(A. 1-10)  we  get 


C  (v   -  V  ) 
-   o 


it_    C^CVy^ 


z  F 


RTC    ( z  ^  -  z  ) 
o     +     - 


(A. 1-11) 


Since  the  solution  is  electrically  neutral 


z^C^  +  z_C_  =  0 


z  C 


Hence 


z  C 


'  + 


+ 


z^   -    z. 


z  C 


-  z 


C_^  +  C_ 


V  C 


v^C  +  v_C 
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V   +  V 


V 

-^  (A. 1-12) 


Substituting  Equation  (A. 1-12)  into  Equation  (A. 1-11),  we 
get 

it_    C  CVy  v_ 

C_(v^  -  v_)  =  -  ^  +   ^,^TC    ^      (^-1-1^) 

o 

Hence 

^JC  it_ 

C_v_  =  -  ^^CV,^  +  ^+  C_v^      (3.3-6) 
o 


Similarly,  it  can  be  shown  that 


APPENDIX  2 
DERIVATION  OF  EQUATION  (3.3-10) 

The  electrochemical  potential  of  the  electrolyte  is 
given  by 


vRT!in(Cf^a°) 


(A. 2-1) 


where  f^  is  the  mean  molar  activity  coefficient  of  the 
electrolyte  and  a^  is  a  proportionality  constant.   The 
gradient  of  the  electrochemical  potential  of  the  electrolyte 
is  then  give  by 


Therefore 


Thus 


Vy, 


vRT  |—  anCCf^a^) 
9x       ±  ±' 


VRT 


Cf  a 

+  + 


vRT 
Cf  , 


3C  ^  ^  0  ^^: 


9f 


f±  ^  ^  8C^ 


CVy 
_j 

vRT 


1  + 


1  + 


VRTC, 


CVu 


C9f 


f^9C 


92.nf^ 
9TnC~ 


c. 


vc 


vc 


vc 


1  + 


dP.nf^ 
d«,nC 


(A. 2-2) 


VC 


(A. 2-3) 
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The  measured  diffusion  coefficient  of  the  electrolyte,  D, 
is  related  to  V   by 


d£nY  + 


D  =  V{1    +    vM  m) (1  +  ,„ 

o        diinm 


(A. 2-4) 


V  =  V    +  V 


m  = 


P  - 


j  =  0 


C  .M. 


C  M 
o  o 


Therefore 


D    =    V 


=    V 


1  +  (v^  +  v_ 


C   +  C^  +  C 
o     + 


M 


O  C  M 
o  o 


(A. 2-5) 


1  + 


dJiny^ 
dUnm 


C 


1  + 


df.ny^ 
dinm 


=  V 


1  + 


dSlny^ 
dJinra 


(A. 2-6) 


Hence  we  get 


D 


1  + 


dJ,nf 
dJinC 


1  + 


d£nY_^ 
d2,nm 


VC 


(A. 2-7) 


The  measured  diffusion  coefficient  of  the  electrolyte,  D, 
is  also  related  to  V   by  the  following  relation 


D  =  — 


1  +  vM  m 
o 

1  +  Mm 

_ 

1  - 


1  + 

dg,np 
d£nC 


de.nf^ 
dInC" 


(A. 2-8) 
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From  Equations  (A. 2-4)  and  (A. 2-8),  we  get 


1  + 

d£nf^ 

de.nC 

1  + 

d£nY  + 

(1  ,  ,^)   a  -  11^) 


d£nm 


= 

1  +  ^'^ 

CM 
*■      o  o 

J 

.- 

d£np 
d£nC 

= 

1  +   MC   _ 

CM 
o  o 

[  1  +  ^M 
^      Po  J 

d£np 
d£nC 

= 

1  4-   ^C   - 
CM 
o  o 

pd£np 
p  d£nC 

= 

1  - 

pdlinp     MC 
I  Po^^nC    p^  , 

p  = 

C  M 
o  o 

+  MC 

d£np 

( 

L 

g'c  (^o^o  ^  ^^^ 

dC 

:  M 

+  MC 

(A. 2-9) 


o  o 


3C 


M 


p   o  9C 


o  _!_  M 


Multiplying  each  side  of  the  above  equation  by  C ,  we  get 


8C 


d£np 
d£nC 

=  ^  M     ° 
p   o  9C 

_^  MC 
P 

c^pac 

MC 
P 

p  PnC 
p£nC 

MC 
P 

(A. 2-10) 


Substituting  Equation  (A. 2-10)  into  Equation  (A. 2-9),  we  get 
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1 

+ 

d^nf^ 

d£nC 

1 

+ 

dK,nY  + 

1  - 


di?.nm 


1  - 


P„ 


d^nC 


p  dSinC 
o o    MC 

p  d£nC     p 


d2.nC 


o  _^  MC 


MC 
P„ 


MC 
P 


o  J 


=  1  - 


dlinC 
c 

d^nC 


(A. 2-11) 


Substituting  Equation  (A. 2-11)  into  Equation  (A. 2-7),  we 
thus  get 


^JWTC-   ^"^e  =  ° 


1  - 


d2.nC 
c 

d£nC 


VC 


(3.3-10) 


APPENDIX  3 

CONVERSION  OF  D^  TO  D* 

(VOLUME  AVERAGE  FRAME  DIFFUSIVITY  TO 

MOLAR  AVERAGE  FRAME  DIFFUSIVITY) 


Since  the  diffusion  data  that  we  have  on  KOH  from 
Bhatia  are  in  the  volume  average  frame  only,  they  have  to 
be  transfonred  into  the  molar  average  frame  to  be  used  in 
our  equations.   The  diffusion  flux  in  the  volume  average 
frame  can  be  related  to  that  in  the  molar  average  frame  in 
the  following  manner.   For  a  binary  system,  from  the  defini- 
tion of  volume  average  frame, 


V  V 

V  A  ,        B 

^  =  Pa^a  m^  ^  Pb^b  M^ 


(A. 3-1) 


•V       ,        V, 
^A  =    Pa^^A  -  ^  ) 


(A. 3-2) 


Substitution  of  Equation  (A. 3-1)  into  Equation  (A. 3-2)  gives 


.V 

^A  =   Pa 


^A  - 


Pa^a  m-  "■   Pb^b  m- 

A  B 


=  Pa  t^A  -  (^A^A^A  ^  ^V^B^B^  ^ 


=  Pa^a  (^  -  ^a^a)  -  Pa^b^b^b 


(A. 3-3) 


Dividing  Equation  (A. 3-3)  by  M-  gives 
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^1=    Va^^    -    ^A^A^     -    ^aSVb  (^•^-^) 


Since      N^  =    C^v^  =    J*    +    C^v*    , 


^B    =    S-B    ^    ^B    ^    S^*     ' 


^"^       ^A^A    ^    S^B    =    1     ' 


Equation     (A. 3-4)    becomes 


^l    =     (^A    ^    ^A^*)     ^B^B    -     (Jb    ^    ^B^*)     C^^B 
-     ^Vl    ^    ^A^B^*    -    ^A^B    -    Vb^*^     ^B 

=    ^Va  -  Vb^    ^b  (^•^-^) 


since  J*  +  J*  =  0,  Equation  (A. 3-5)  becomes 

A      rJ 


^1=^  ^A  (^A^  S^  ^B  (^•^-^) 


Since  J^  =  -d"^VC   and  J   =  -D  VC^,  (A.  3-6)  becomes 
-d"^VC-  =  -D*  iC^+C^)    V„VC^  .   Therefore,  D*  =  


(VS^^B 


APPENDIX  4 


CALCULATION  OF  C2(r,x)  FOR  A  GIVEN  x 


The  flux  equation  for  C-,  is  [26] 


d  C-    D^  dC^    ,   di 

D   I   +_J._^_i ^  =    0 

3,2     r  dr     RF  dx 
dr 


(A. 4-1) 


Substitution  of  Equation  (5.1-1)  into  Equation  (A. 4-1)  gives 


d^C     D   dC 

D-3  T^    +    —   ^    -    C, 

3,2     r  dr      3 
dr 


ai 


nFC. 


C^(x) 


,o 


exp 


RT 


ai 


nF 


C^(x) 


exp 


RT 


=  0 


(A. 4-2) 


Let 


ai 


C^ 
o 


nFC°D2   C. 


exp 


F(*3-^,) 


RT 


(A. 4-3) 


and 


C 
^3 


C3    C3 


C^(x)  C^ 
l_C°  C^(x) 


exp 


r2F(*s-*£) 


RT 


(A. 4-4) 


Then  Equation  (A. 4-2)  becomes 


d  C      dC' 

3    13     2 

^    +    -  -^   -  q^C^  =  0 


dr 


r  dr 


(A. 4-5) 
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The  following  boundary  conditions  can  be  used  in 
solving  Equation  (A. 4-5), 


B.C.I 


At  r=r  ,  Cl 
o         3 


=  c";  -  C. 


C,  (x)   C 


C^Cx) 


exp 


2F(*^-*^) 


RT 


C 
o 


where  C,  is  the  equilibrized  concentration  at  the 
gas-electrolyte  interface. 


dc: 

B.C. 2:   At  r=r  +6,   :3— ^ 
o      dr 


=  0 


The  solution  to  Equation  (A. 4-5)  is 


C-  =    A^I^(qr)  +  A^K^Cqr) 


(A. 4-6) 


where  A,  and  A^  are  constants  and  I  (qr)  and  K  (qr)  are 
1       2  o  ^         o  ^ 

Bessel  functions  of  zero  order.   Applying  B.C.'s,  we  have 


at  r=r  ,    C  =  B, I  (qr  )  +  B^K  (qr  )  (A. 4-7; 

o      o     lo^o      2oo 


at  r=r  +6,   0  =  B  I  (q(r  +6))  +  B^K, (q(r  +6))   (A. 4-8) 


11 


2''1 


From  Equations  (A. 4-7)  and  (A. 4-8)  we  have 


A-,  = 


c;K^(q(r^+6) 


1    K^(q(r^+5))I^(qr^)  +  K^ (qr^) I^ (q (r^+6 ) ) 


and 


A^  = 


CJ^^{q{r^+6)) 


K^(q(r^+6))I^(qr^)  +  K^ (qr^) I^ (q (r^+6 ) ) 
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Substituting  these  into  Equation  (A. 4-6)  gives 


c;  = 


c' 

o 


3    K^(q(r^+5))I^(qr^)  +  K^  (qr^)  I^  (q  (r^+6  )  )  "" 


[K^(q(r^+6)  )I^(qr)  +  I^  (q (r^  +  6 )  ) K^ (qr )  ]  (A. 4-9) 


or 


C.  -  C 


o  ^1^^)  ^4 
^  C°  C^(x) 


exp 


RT 


e   o  ^1^^^  ^4 


exp 


2F{(t)g-<J)j^) 
RT 


DEN 


X 


[K^  (q(r  +6))I^(qr)  +  I,  (q (r^+6 ) ) K  (qr) ]    (A. 4-10) 


where  DEN  =  K^ (q (r^+6 ) ) I^  (qr^)  +  K^ (qr^) I^ (q (r^+5 ) ) 
Dividing  Equation  (A. 4-10)  by  C°  gives 


'A 


C°  C^lx) 


^^P  RT ^ 


,o 


C^(x)  C° 


exp 


2F(^^-*^) 


C3    C^  C4(x) 
[K^(q(r^+5))I^(qr)  +  I-^  (q  (r^+6  )  )  K^  (qr)  ] 


DEN 


RT 


(A. 4-11) 


Takinq  the  average  using 


C3(x)  = 


f    2T:rC-,  (r  ,x) 


•rrr 


dr 


gives 


^o 


C^(x)  C^ 
C°  C^(x) 


exp 


2F(.}.g-<})(,) 


RT 


2r 


o 


C  (x)  C^ 
1  -  X 


C^  C^ix) 


exp 


2F(*^-^^) 


RT 


q (DEN) 
[I^(q(r^+6))K^(qr^)-K3^(q(r^+6))I^(qr^)] 


2r  6  +  6' 
o 


(4.5-2) 
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C3(x) 
With  this  average  value  of  —^ ,  Equation  (5.1-1)  becomes 


di 
s 

dx 

=   ai 
o 

C3(x) 

exp 

'^1 

RT 

C-^(x) 
-      p           exp 
^1 

rF(4)3-(|)^) 

RT 

(4.5-3) 

APPENDIX  5 

COMPUTER  PROGRAMMING  FOR  THE  COMPUTATION  OF 
CURRENT  DENSITY  DISTRIBUTION 


The  conputer  programming  for  the  computation  of  current 
density  distribution  is  listed  below: 
C 

t: 

r 

i"     M/i  ri-ii  MATir  A[    "iinriiiKiQ  of   a   fufl  hf.i  l    FiFCTRnne 
c 

dp'Fmsti'm     rKiiHfpin,    rMp(i(?n),    ph\(.?o),    p^s(2c\) 


Rr  61       t^ J  ,    t:?,    K^, .  Ka 


L' 
I" 

I"       VAPlATTnH    CF     T'-'F      /il  THOpi'' iMf:    CIPi^EmT 
C 

,IM     =      )? 

DM        ^(M)        J     =     ]  ,     JM 

ClIRAPp     =     f'.jon     A      T 

rURf A1  f  n      =     0,0 

F»-'|    A'^'^     :     0  .SO     +     0  .U'^S      A     J 

ni  L  rp^'    =    n.dS 
L  V    =    I  0 
111      Hfi      ado      I     =    1  ,    I  ►' 

P  ML  ASS     =     t'HI   tSS     -     Dl  L  TPH 
PMI    f  1  1     =     PHI    aSS 

FHsn  )   =   o.y 

ni  PHS!    =  pHp  (  n   ~  PHI.  f  1 ) 

c     Mi'MFi?ir./i    T'rrr^.PMTriM 
c 

N    =     'w' 

THICK      =     0,0] 

w  I  n  r  H   =   T  H  T  r  K   /  k 
c 

A  (<  n  1    =  5  a  ,  H  ^    A    n  I  p  H  s  [ 

A  D  G  ?    =  A  n  (*  1     A    ?  .  0 

p  M  s  C  T  1  =     iv  T  r^  T  H     i     ] 

AA     =     f  CMIRAPP     -     CHRfAT(T)     )     A     ^S9.0 

CNKPH    =    r.^o^Ci) 

I' 

I'        PAPFTAI      Mill   A I      vniliMf-     (IF     l-AT^R 

l: 


pARV'i   =   (i.uiei'ro?   +   rjKfiH   a    (   7,37a6"i*i  o.  a*  f-^,)    _ 

C  ^^'K(1H  X     (     a.7^"5l  ^A  1  0  .  A  A  (-r-,;)  _  CMKHH  a  (  a, 67^09 

r  A10,>AC-61  -  1  .Ab^AaA 1 n.  A  a (-7)  A  CNKdH  ))) 


C       CHMCFfi  TPM  THK,    HF     W  A  T  F  R 


c 


CHpiMTI    =    '^'4.75S^    -   (:^'KnH    a    (    a.5in'^    +   CMKfn-    a    f    "^.1709* 
r,  lo.AAf-.M    -    1  .^^i3m  0.  AA  (-n    A    cNKrM^    )) 

CMH?n   -    CH?(i(I ) 
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i: 


PftUTl  Al      r'lil  *l      vn[  ij-^T    (IP     I-  'Hi 

PAkV  =  f  1.00  -  rf'HPn  *  PARVU)  /  CMKOH 

Di.nvnv    =    tA    *    (    i.ii    -   r.fjKriHxPARv    )    *    (    i.n    -   ^^!K^H    x 

r  r     PAI7V     -     PiRVU     )1      /     (      1  .  fl     -     CMKCiH     A     PAPV     > 

r.    =    f,rf)r    f    i.7ni,'i:<,    a    rNH?n   a   rvnc-ArK;])    5 


Rt    =    r: 
'r,>    =    n 

CAl  1 
fAl  I 


*  !  .  0(1      A      1  (I  .  A  «•  f-r^j  5 

*  1  .  O^J      A      1  II  .  i  *      (-b) 


rfont    c    RSI  T 1 1  ,    RSI  xi  1  ,    ri  ,    ^    5 

FFO/IA     (     RSITI?,     HSlKl?,     R^J,     ?     ) 


DM)    =    F;.-Ln?     A    n.siKH     -    bSLKl?    A    RSLIll 

r.rM   =   p SI. KIP   A   H.si  III   +   MSI  Ko    A   R:uiir 

CdVYCT)     =     ol'JOOiqPl     *     (     (     CMKOH/CMH?!!     )     *     pvpfARr;?^ 

r  f    P.^iSi  A 1  (I .  A  A(61    A    f    1.0    -    S.Pa«    A    f    r;^"<n'^ 

C  /  r  ^,  M  ^-^  fl    ^     ^    1^  X  P  (  A  R  n  p  1     ^     A    D  ^J  1 1    /     (    f,-    a    n  E  ^'    )     ) 

rMovY  =  cdvYtn 

CmTilT  n  V  T  rv    Of    FL  F  CTRfil,  Yir 
Cf  INT  I  If 


c 


Di  ir 


r. 
r 


-0„S^P7C-,     +    C'^Krih     A     c     0.60H7S7     -    C'lKflH     a 
C    o.Ob'^apnh    -    r,  .0(}??9'hbn    a    rMKiiH    5    ) 

rrlKi'lri     A     r      16P1.3^J57«     A     CNKHH     -     f^SP^.P^TM     )     / 
(    Pl^OJij6U'5    -    O.xf.H    A     (    f^h?.b,2P>7n    -    Pin. 67187 
A    C^KPH    1     ) 


Dt-LPHi    =  f    f  cuRCATcn   -  njRAPp   )   /  ciifiniir  -   r   o.o'^tui 

C  /    O-'^HM    1    A    r    o.<^Pa    +    P.^1    A    CMKnH/rNnpn    )    a 

C  nrnvny    A  '  DI  A  ni.f     )     A    winiH 

pui  f  T  Ml  =     '^t't  n)     +     OPLPhl. 

A|  pt'i     =     /4  ~Si)     A     ]  0  .  A  A  (  +  S  ) 
XLUfRn    r    ^.i)     A     10.AA(-fll 

oircrrCT"!    =    alpha    a    yrHct'O   a    c    o.7aP'i    a   p'^inxY    a    r^iHPri 

r,  A      Kync_ARnn      -      C      CNKOH/^,  ,7979      )      A 

r  e;  X  p  f  A  R  r,  n    )    A  v.  T  n  T  H 

niprATCnn   =  cmrp^tct)    +•  ni.cf'CTci) 

RHP  =   1 1100 

[u'l  PF-s   -    f-ru'-r  AT  f  i  +  n    /   fnm   )    a   width 

p  H  s  f  T  >  n    =    P  ^'  "^  c  T )    +    n  F  I  P  h  ? 

ni  PHSI      =    PMSfI  +  !l     -    P'-'LCI  +  I) 

CAI  r.PI  ATld'^i    (IF    KHM    rf1Nf:r>:TRA  rifU;    hy    tmf 

niiAf,RArTr   pi'i-^r -^^i'tt.'    MFTurn- 


r 
c 


tu:   =   I.. iriH    A 
K) 


A  A 


'^H  A  f  I.n  -  C^JKdHAp  ARV  )  A  f  1.0  -  CNt-'dH  A  ( 
PARV  -  P.*RV/1  )  )  /  (  1  .  (j  -  CMK(lh'*(  PARV  -  ?.n 
PARV!   1  ) 


I) 


HH     A     (      1.0     -     C     CMKilh     +      <l/P,(l     "1      *     RARV     )      A      (     1,>, 
-     r     CMWfiH     +     Kl/^.O     )      A      f     PApV     -     PARV/l     ))/(].() 

(    vy.^ny    +    k  1/2,0    )    a    (    i'APv    -   2,c,    a    pa,pv«    )    ) 
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lOU 


11 


a  0  0 

cfOO 
101 


lOi? 


lO'j 

iou 


Ki     =     H«     *     f     1.0     -     (     C^'KnH     +     KP/?.()     I      *     PAPV     )     *      (     1  .  (i     - 

r.  f   rnKni-'   +   K?/?.o    )    A    (   pARv   -   PAru"!    i    )    /    r    i.n   - 

r  (     CNKOH     +     Kc^/P.O      )      A      (     pA';V     -     2.0      A     DARV/i      1      5 

K/1    =    f^P     A     (     1.0     -     (     C^Knl■>    t    K5     1     A     PAI^v     1     A     (     1.0     - 
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